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A.  Abstract 

The  utility  of  the  plasma  focus  as  an  opening  switch  was  analyzed  on  the 
basis  of  our  previous  experimental  observations  and  those  on  other  Mather-type 
devices,  A  load  resistor  has  recently  been  included  in  the  opening  switch  cir¬ 
cuit  model  and  the  analysis  indicates  that  the  typical  pinch  resistances 
measured  (~0.1ft)  are  too  small  to  drive  typical  load  resistances  encountered  in 
pulsed  power  in  terms  of  %  energy  transferred  (<2%)  and  the  time  required 

for  the  transfer  (>200ns).  A  scaling  law  for  pinch  stability  has  been  derived, 
which  is  consistent  with  our  data  and  that  from  several  other  experimental  devi¬ 
ces;  this  model  predicts  that  a  plasma  focus  could  be  designed  to  permit  longer 
energy  transfer  times  by  increasing  the  anode  radius. 

As  to  the  progress  of  the  two  students  supported  by  the  contract,  Mr.  Frank 
Venneri,  has  successfully  built  and  operated  a  curved-crystal  soft-Xray  sepctro- 
meter  and  he  is  in  the  process  of  writing  his  Ph.D.  thesis.  Mr.  Kevin  Boulais 
completed  the  requirements  for  the  Masters  of  Science  degree  in  electrical  engi¬ 
neering.  Mr.  Boulais'  thesis  on  laser  interferometry  of  the  pinch  phase  along 
with  Mr.  Venneri's  spectrometer  finally  gave  us  a  set  of  state-of-the-art 
diagnostics  with  which  to  study  the  plasma  physics  of  the  plasma  focus.  The 
results  of  these  studies  are  being  prepared  for  publication  . 

B.  Research  Objectives 

As  stated  in  the  proposal,  the  overall  goal  of  the  research  was  to  complete 
the  testing  of  the  plasma  focus  as  an  opening  switch  (PFOS)  in  the  Filippov 
geometry  with  sufficient  time  allotted  so  the  two  graduate  students  could 
complete  the  experiments  required  for  thier  thesis.  Since  the  utility  of  the 
PFOS  concept  has  always  been  uncertain  and  since  theses  must  contain  successful 
results  for  publication  and  to  further  careers,  the  testing  of  the  PFOS  concept 
was  not  chosen  for  either  thesis. 


C.  Status 

While  I  was  able  to  perform  considerable  analysis  as  to  the  merit  of  the 
Filippov  geometry  in  the  opening  switch  configuration  of  the  load  coupled 
directly  to  the  anode  (Fig.  1),  the  experiments  required  to  support  the  student 
theses  tied  up  the  plasma  focus  device  for  the  entire  period  of  funding.  While 
there  is  a  small  chance  the  Filippov  experiment  will  be  performed  this  spring 
while  Mr.  Venneri  is  writing  up  his  Ph.D.  thesis,  there  are  no  new  experimental 
results  at  this  time.  However,  the  experiments  involving  the  student  theses 
were  enormously  successful,  being  sort  of  a  rite  of  passage  for  both  students 
and  launching  them  both  on  promising  careers  as  experiment  research  sceintists. 

In  the  first  section  Cl  the  analysis  and  prospects  for  the  plasma  focus  as 
an  opening  switch  will  be  presented  briefly,  with  the  details  supporting  the 
conclusions  included  as  Appendix  A.  Then  the  highlights  of  the  Francesco 
Venneri' s  and  Kevin  Boulais'  theses  are  presented  in  section  C2  and  C3;  Mr. 
Boulais'  Masters's  thesis  is  included  as  Appendix  B.  A  copy  of  Mr.  Venneri 's 
thesis  will  be  sent  to  AFOSR  as  soon  as  he  has  completed  the  requirements  for 
the  Ph.D. 
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Cl .  Analysis  of  the  PFOS  Concept 

An  equivalent  circuit  for  the  plasma  focus  load  circuit  is  shown  in  Fig.  1. 
The  energy  is  stored  in  the  external  circuit  inductance  L$  which  drives  the 
plasma  to  a  pinch  at  t=0.  The  resistance  of  the  pinched  plamsa  is  represented 
by  R0(t)  taken  to  be  a  step  function  beginning  at  t=0  when  S2  is  closed.  For 
the  fast  transients  considered  here,  it  is  assumed  the  capacitor  bank  has  a 
negligible  effect. 

Previous  solutions*  included  only  a  load  inductor,  but  recently  the  solutior 
was  extended  to  include  a  series  load  resistor,  RL.  This  permits  the  model  to 
be  applied  to  more  realistic  loads  including  transmission  lines,  or  pulse 
forming  networks  which  a~e  commonly  used2,3  to  match  Marx  generators  with  beam 
diodes.  The  impedance  of  these  networks  are  typically  in. 2 >3 

The  analytic  solution  is  as  follows.  The  load  current,  II  is  given  by: 
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From  this  the  voltage  drop  and  energy  transferred  to  the  load  resistor,  V|_  and 
E|_  respectively,  can  be  calculated  and  the  time,  t^,  of  the  peak  power  transfer 
to  the  load  can  be  found. 

From  our  previous  experiments4,  where  the  peak  load  current  I|_  ~  14kA,  Ls  = 
17nH,  L0  =  29.4nH,  Ll  =  30nH,  and  R|_  =  0.8ft,  one  finds  that  the  solution  best 
fits  the  data  if  R0  ~  0 .07 sis .  This  value  is  in  the  range  observed  in  other 
plasma  focus  experiments.5 

If  one  assumes  Ro  ~  0.1ft  independent  of  the  device  and  varies  Rl  for  the 
Limeil5  Filippov  plasma  focus  (which  is  of  more  optimal  design)  where  Ls  =  17nH, 
L0  =  21.9nH,  Ll  =  6nH  one  can  find  the  peak  II,  Vl,  ElUm),  Ej_ (cc)  and  tM  for 
each  Rl  from  equation  1.  This  is  shown  in  Fig.  2. 

As  can  be  seen,  the  ElUm),  ElC")  and  iLpk  all  suffer  a  decrease  with 
increasing  Rl,  especially  for  Rl  >  0.1ft.  While  the  peak  voltage  increases  with 
Rl  as  desired,  these  effects  are  overcome  by  the  decreases  in  ILpk  and  El(°°). 

The  latter  being  less  than  2%  of  the  total  inductive  energy  available  E0  (E0  = 
1/2(Ls  +  L0)I02  =  63kJ). 

The  problem  is  that  the  opening  resistance  of  the  plasma  focus  appears  to  be 
too  low  to  permit  the  transfer  of  a  significant  portion  of  the  initial  energy  on 
a  short  time  scale.  This  can  be  seen  in  the  following.  If  one  assumes  a  fixed 
Rl  of  1.0ft  in  the  previous  example  and  varies  R0  between  0.1  and  10  ohms  one 
finds  the  ratio  of  El('°)/E0  increases  from  0.018  to  0.276  and  t^  decreases  from 
52ns  to  6.3ns  (Fig.  3).  The  fraction  of  the  initial  energy  E0  stored  in  the 
inductances  Ls  and  L0  (E0  =  63k J  for  this  example)  transferred  to  the  load  (by  t 
=  4tM)  is  ~  0.01Eo  for  R0  =  0.1ft  to  0.23Eo  for  R0  =  10ft  with  EL(4tM) >0.1Eo  for 
4tM  <  92ns  for  Ro>1.0ft.  The  rest  of  the  energy  E0  is  consumed  in  R0  during  this 
period. 

Thus  the  plasma  focus  would  be  a  potentially  useful  puisea  power  device  for 
Ro  >  1.0ft  but  unfortunately  our  measurements  are  consistent  with  a  value  of  less 
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than  10%  of  this.  While  at  best  (for  R0  =  10O  over  3/4  of  the  inital  energy  is 
consumed  in  the  switch)  the  fast  energy  transfer  time  (4t^  ~  25ns)  would 
make  it  very  attractive,  for  the  R0  observed  (~0.1ft),  the  transfer  times  are  too 
long  (4tfo  ~  200ns)  and  the  fraction  transferred  too  small  (EL(4tM)~0.01  E0)  to 
seriously  consider  the  focus  as  a  candidate  for  an  opening  switch. 

In  all  this  analysis,  it  was  assumed  that  R0  was  constant  for  several  times 
t|vf ,  however  in  reality  the  enhanced  resistance  only  occurs  over  roughly  the 
duration  of  the  pinch  (or  the  period  of  stability)  t.  For  Mather-type  plasma 
focus  devices  this  period  can  vary  from  10-20ns  for  a  lkJ  device  to  ~150ns  for  a 
device  of  a  few  hundred  kilojoules.  Through  considerable  analysis  of  our 
experimental  data  and  that  from  other  devices,  we  have  developed  a  scaling  law 
to  predict  this  period.  This  analysis  is  useful  to  determine  whether  the  energy 
transfer  could  actually  occur  or  be  interrupted  prematurely  if  T<4tM)#  jhus  -ft 
would  be  essential  in  developing  the  plasma  focus  as  an  opening  switch  espe¬ 
cially  if  Ro>1.0n  and  is  included  here  as  Appendix  A. 

Basically  in  the  analysis  it  is  found  that  the  stability  period  t  =  a  r0/Vr 
where  r0  is  the  pinch  radius,  Vr  is  the  radial  collapse  velocity  and  a  is  a 
constant  of  proportionality;  this  is  consistent  with  experiment.!  Since  empiri¬ 
cally  r0  a  r/\!"3  where  r/\  is  the  anode  radius.  Longer  t  and  hence  transfer 
periods  may  be  obtainable  by  increasing  r0. 

C.2  Francesco  Venneri's  Ph.D.  Thesis  Project 


In  order  to  measure  the  electron  temperature  during  the  pinch  to  determine 
the  Spitzer  resistivity  and  the  speed  of  sound,  a  dispersive  element  is 
necessary  to  avoid  errors  due  to  contamination  of  the  spectrum  due  to  line 
radiation  caused  by  impurities.  For  example,  using  a  filter  technique  and 
assuming  a  free-free  bremsstrahlung-type  spectrum,  plasma  focus  researchers  have 
tended  to  estimate  the  electron  temperatures  to  be  in  the  multi  keV  range 
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whereas  using  laser  scattering  they  tend  to  be  in  the  sub-keV  range.?  In  our 
device  using  filters  and  assuming  a  bremsstrahlung  spectrum,  we  have  measured 
electron  temperatures  to  be  in  the  300-500eV8  range  but  the  presence  of  line 
radiation  can  not  be  ruled  out. 

A  dispersion  element,  such  as  a  crystal  spectrometer,  permits  the  dispersion 
of  the  radiation  as  a  function  of  a  wavelength  so  the  spectrum  can  be  observed 
photographically  and  the  regions  of  line  radiation  and  bremsstrahlung  can  be 
identified  and  separated.  Thus  the  radiation  can  be  analyzed  using  the  proper 
model  for  that  region  of  the  spectrum  and  an  accurate  value  of  the  electron  tem¬ 
perature  determined. 

For  Mr.  Venneri's  thesis  project,  it  was  decided  to  construct  a  state-of- 
the-art  curved-crystal  spectrometer^  for  this  purpose.  In  the  following  sub¬ 
sections  the  spectrometer  will  be  described  and  some  typical  experimental 
results  presented.  How  these  results  have  been  applied  to  the  project  is 
illustrated  in  Aooendix  A. 

A.  Curved  Crystal  Spectrometer 

A  crystal  spectrometer  uses  the  diffractive  properties  of  crystals  to  iso¬ 
late  wavelength  in  a  distribution  of  x-rays,  so  that  they  can  be  detected  indi¬ 
vidually.  Due  to  the  regularities  of  crystal  lattices,  x-rays  of  the  same 
wavelength,  scattered  from  different  crystalline  planes,  experience  constructive 
or  destructive  interference,  depending  on  the  differeneces  in  their  path 
lengths. 

The  interference  is  such  that,  for  crystals  with  a  sufficiently  large  number 
of  parallel  planes,  only  the  x-rays  with  wavelengths  obeying  the  simple  Braggs' 
rule  are  constructively  diffracted  (reflected)  at  a  certain  angle  (x)  (Fig.  4): 


Where:  n  =  Order  of  the  diffracted  beam  (equal  to  the  number 

of  the  first  plane  under  the  surface  from  which 
the  wavelength  is  constructively  diffracted), 
d  =  interplanar  spacing  of  the  crystal. 

1  =  difracted  x-ray  wavelength, 
x  =  angle  of  diffraction. 

When  analyzing  a  pi  used  source  such  as  a  pinched  plasma,  a  plane  crystal  can 
offer  only  a  limited  range  of  diffraction  angles  to  an  incoming  (diverging)  x- 
ray  beam,  and  no  range  at  all  if  the  beam  is  collimated  parallel.  To  detect  a 
wider  portion  of  the  spectrum  a  curved  crystal  should  be  employed. 9 

In  the  x-ray  spectrometer  developed  for  the  Dense  Plasma  Focus,  the  curved 
crystals  are  elliptically  bent  with  the  plasma  filament  (pinch)  at  the  far  focus 
of  the  ellipse.  The  cuvature  in  the  crystal  produces  a  continuous  range  of  » 
incidence  angles  between  the  incoming  beam  and  the  crystal  planes.  Different 
frequencies  are  constructively  diffracted  at  different  angles,  then  converge  to 
a  common  crossover  point  (the  near  focus  of  the  ellipse)  and  finally  diverge 
towards  the  detection  circle. 

The  beam  is  unrestricted  (diverging)  in  the  horizontal  plane,  but  collimated 
in  the  vertical  direction  (pinch  axis)  to  exclude  viewing  of  the  electrode 
region.  At  the  crossover  point,  a  collimating  vertical  slit  (parallel  to  the 
pinch)  is  placed  to  reduce  the  background  radiation  from  regions  outside  the  far 
focal  axis  of  the  ellipse  (location  of  the  pinch). 

On  the  detection  circle,  an  array  of  solid  state  detectors  (PIN  diodes)  pro¬ 
vides  time  resolved,  space  averaged  readings  of  the  X-ray  emission  at  selected 
spectrum  locations. 

An  aluminized  Kinfol  vacuum  window  (O.lu  aluminum  deposited  on  2u  Kimfol) 
filters  visible  and  UV  light  out  of  the  spectrometer,  which  is  kept  under  oil- 
free  vacuum  to  prevent  attentuation  of  the  x-rays,  contamination  of  the  detec¬ 
tors  and  breakdown  of  their  bias. 


Different  crystals  can  be  used  in  the  spectrometer,  each  with  its  own  dif¬ 
ferent  interplanar  spacing  asnd  useful  frequency  range,  not  to  mention  bending 
and  mounting  techniques.  The  geometric  resolution  of  the  spectrometer  depends 
on  the  match  of  detector  and  crystal  geometry.  At  a  detection  circle  of  15  cm 
in  radius  the  resolution  is  between  .5  and  2  ev/mm,  depending  on  the  eccentri¬ 
city  of  the  crystal  (dE/E  =  1/1000  for  KAP  and  PET,  1/20000  for  Li F ) . 

B.  Temperature  Measurements 

The  electron  temperature  of  our  DPF  was  determined  as  a  function  of  time, 
using  the  curved  crystal  spectrometer  (Fig.  5). 

After  a  line-free  region  of  the  spectrum  was  located  (energy  regions  above 
the  highest  line),  the  bremsstrahl ung  emission  was  monitored  with  three  PIN 
diode  detectors.  This  allows  to  determine  the  temperature  as  a  function  of  time 
in  the  pinch. 

The  temperature  is  in  fact  proportional  to  the  logarithmic  slope  of  the 
bremsstrahl ung  spectrum  and  the  use  of  three  detectors  allows  the  verification 
of  the  effective  thermal  nature  (as  opposed  to  beam  target)  of  the  monitored 
continuum. 

This  method  is  applicable  (with  high  efficiency  crystal  geometries)  even  to 
pure  hydrogen  or  deuterium  discharges.  There  is  no  need  to  use  a  seed  gas  to 
enhance  the  x-ray  emission,  unless  it  is  of  course  desired  to  investigate  the 
spectrum  of  such  seed  gas. 

For  hydrogen  and  deuterium  discharges  the  temperature  was  determined  using  a 
100  mm  x  20  mm  PET  crystal  of  eccentricity  0.806  placed  at  a  focal  length  of  65 
cm  from  the  Plasma  Focus  pinch.  The  energy  bands  selected  for  this  application 
were  1468  eV  (+-5  eV)  and  1670  eV  (+-6  eV).  The  spectrometer  was  independently 
calibrated  using  a  foil  filter  technique  (filtered  detectors  mounted  at  the  end 
of  a  4  m  collimating  tube)  and  a  second  crystal  spectrometer  (flat  PET  crystals 
at  1592  eV  and  2006  eV).  The  relative  efficiencies  of  the  PiN  diode  detectors 
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were  also  determined  employing  the  long  collimating  tube  to  view  the  same  region 
of  the  plasma  pinch  during  a  series  of  calibration  discharges. 

The  measure  temperatures  range  between  200  eV  and  400  eV  (Fig.  6). 

The  temperature  obtained  using  this  technique  is  very  nearly  the  maximum 
temperature  obtained  at  the  most  dense  part  of  the  pinch,  as  it  develops  as  a 
function  of  time.  In  fact  the  dependence  of  the  bremsstrahlung  radiation  on  the 
square  of  the  density  and  exponentially  on  the  temperature  heavily  weighs  the 
hottest  and  most  dense  regions  of  the  pinch  as  the  radiation  is  averaged  over 
the  viewed  emitting  region. 

These  results  are  the  first  using  this  type  of  spectrometer  on  a  plasma 
focus  and  generated  considerable  interest  during  two  presentations  at  plasma 
physics  meetings.  10-H 

C.3.  Kevin  Boulais'  Masters  Thesis 


Although  we  have  been  operating  our  plasma  focus  since  1978,  we  had  no 
measurement  of  the  plasma  density  or  its  profile  until  this  year  when  Mr. 
Boulais'  laser  interferometer  came  on  line.  With  such  an  instrument,  one  can 
obtain  a  snapshopt  of  the  density  profile  in  the  plasma  focus  (see  both 
appendices)  at  one  instant  (At<5ns)  per  shot.  Through  n  series  of  runs,  we  have 
determined  that  the  dynamics  of  the  plasma  (as  opposed  to  the  neutron  yield)  are 
reproducible  enough,  so  that  a  reasonable  understanding  of  the  events  taking 
place  can  be  obtained.  This  is  clearly  shown  in  Appendix  A  where  the  decay  in 
the  number  of  particles  per  unit  length  can  be  used  to  obtain  an  estimate  of  the 
magnitude  of  the  axial  flow  velocity.  This  velocity  appears  to  be  great  enough 
to  explain  the  stability  of  the  pinch.  Before  Mr.  Boulais1  interferometer  came 
on  the  scene,  we  could  only  speculate  about  this. 

Mr.  Boulais'  achievement  is  clearly  remarkable  since  no  one  on  campus  had 
any  experience  in  this  area  and  he  only  interacted  with  experts  briefly  at  a  few 
meetings.  Never-the-less  the  interferograms  obtained  are  remarkable  for  their 
clarity. 
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Mr.  Boulais  is  now  pursuing  ,iis  career  investigating  microwave  generation  by 


relativistic  electron  beams  at  the  Naval  Surface  Weapons  Laboratory  in  Silver 


Spring,  Maryland  and  hopes  to  continue  toward  a  Ph.D.  in  electrical  engineering 
at  the  University  of  Maryland  next  fall  (1988). 


For  the  details  of  his  achievement,  his  M.S.  thesis  is  reproduced  as 


Appendix  B. 
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Figure  la. 


The  operation  of  a  plasma  focus-load  configuration  with  a 
load  tap  behind  the  anode.  The  switch  S2  is  closed  at 
pinch  time  trapping  the  magnetic  flux  associated  with 
plasma  focus  inductance  (L0)  in  the  load  circuit.  When 
the  flux  6<f>  subsequently  leaks  out  the  current  II  is 
driven  through  the  load  impedance  Z|_. 


Figure  lb. 


The  circuit  modeled  in  the  calculations  presented  in  this 
section  where  the  plasma  focus  resistance  is  represented 
hy  R0(t ) .  The  load  resistor  R|_  has  been  added  to 
generalize  the  model  so  it  may  be  more  readily  compared 
with  the  experiment. 
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Figure  2.  Effect  of  varying  the  magnitude  of  the  load  resistor,  Rl, 
for  the  Limeil6  Filippov  plasma  focus  where  R0  =  O.in,  Ls 
=  17nH ,  L0  =  21.9nH,  Ll  =  6nH  and  E0  =  63kJ  (Fig.  1).  El 
(«),  Itpk»  tM.  ElUm)  and  VL  are  the  maximum  energy 
transferred  to  Rl,  the  peak  current  in  Rl,  the  time  of  the 
current  peak,  the  energy  transferred  to  Rl  by  tM  (Ro  is  a 
step  function  beginning  at  t=0  when  II  =  0).  VLpk  is  the 
peak  load  voltage. 
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Figure  3. 


Effect  of  varying  R0  on  EL(»),  EL(tM),  Q(tM) ,  Ipk,  Vpk, 
and  for  R|_  =  l.Ofl.  Q(tn)  is  the  reactive  energy  stored 

in  the  load  at  tM.  The  rest  of  the  circuit  parameters  are 
the  same  as  in  Fig.  2. 
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Figure  4.  Illustration  of  Bragg's  dispersion  law  for  a  crystal 

spectrometer.  By  choice  of  different  d  (i.e.  different 
crystal  types)  different  regions  of  the  soft  xray  spectrum 
may  he  investigated. 
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Figure  5. 


if 


Elliptical  curved  crystal  spectrometer  used  to  determine 
the  electron  temperature.  The  plasma  column  is  at  one 
focus  of  the  ellipse  and  the  cross  over  point  at  the 
other.  Pin  diodes  and/or  photographic  film  on  the 
detection  circle  are  used  to  detect  different  spectral 
bands  of  the  soft  x-ray  emission. 
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6.  Typical  time  dependence  of  the  electron  temperature  of  the 
pinched  plasma  obtained  with  the  curve  crystal  spectro¬ 
meter  for  Po  =  3  torr  (deuterium  and  I  =  550kA. 

Comparison  with  simultaneous  laser  interferograms  indicates 
the  first  peak  corresponds  to  first  compression,  the  dip 
to  radial  expansion  and  the  second  peak  to  second 
compression. 
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Abstract 

By  use  of  a  streak  camera  aligned  either  along  the  axis  of  the  pinched 
plasma  or  at  90°  to  it,  the  distance  along  the  axis  between  first  and  second 
compression,  h,  and  the  minimum  radius  of  the  pinch,  r  ,  have  been  measured  as 
a  function  of  the  plasma  current,  I,  and  filling  pressure,  p  .  The  plasma  current 

was  measured  by  magnetic— pickup  loops  calibrated  to  35MHz  with  Helmholz  coils. 

The  results  are  that  h  «  IXpQ^  where  x  =  — 0.02±0.12  and  y  =  0.26±0.13  and 
that  rQ  x  Iap0^  where  a  =  — 0.18±0.11  and  b  =  0.2±0.03.  Thus  within  the 
experimental  uncertainties  h  =  a  rQ  where  a  =  8.6±2.4.  These  results  may  be 
correlated  to  the  stability  period  r  by  use  of  the  zippering  effect;  the  latter 
predicts  the  radial  collapse  velocity  vf  should  be  directly  related  to  the  speed  v 
of  the  compression  wave  traveling  along  the  column.  These  latter  two  quantities 

are  found  to  exhibit  the  same  dependences  with  I  and  pQ  within  less  than  10% 
experimental  uncertainties  in  their  power-law  slopes.  The  empiracle  scaling  of  r 

with  pQ  and  I  found  here  is  consistent  with  results  of  other  plasma  focus  devices 
and  is  consistent  with  ideal  MHD  stability  theory.  Experimentally,  r  is  found  to 
be  enhanced  over  that  of  a  straight  cylinder  by  a  factor  a'  —  2.7±0.6.  Two 
simple  models  consistent  with  ideal  MHD  theory  could  explain  this  enhancement. 
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A  Scaling  Law  for  Macroscopic  Stability 
of  the  Mather— type  Plasma  Focus 
G.  Gerdin,  F.  Venneri,  and  K.  Boulais 

I.)  Introduction 

1  t> 

Stability  of  the  plasma  focus  ’  or  non  cylindrical  Z  pinch  is  of  interest 
for  both  theoretical  and  practical  reasons.  Theoretically  the  stability  of  the  Z 
pinch  represents  a  classic  problem  in  magnetohydrodynamic  (MHD)  theory,  which 
predicts  the  column  should  be  unstable  to  the  MHD  equivalent  of  the 

Q  J 

Rayleigh— Taylor  modes  ’  and  thus  presents  an  opportunity  for  a  direct 
comparision  between  theory  and  experiment.  From  a  practical  point  of  view  Z 
pinches,  both  compressional  and  gas-embedded,  could  represent  a  relatively 
inexpensive  approach  to  thermonuclear  fusion  if  the  period  of  stability  can  be 
made  long  enough.  The  plasma  focus  represents  a  cross  between  these  two 

approaches  to  Z  pinches  being  formed  in  a  similar  manner  to  the  compressional  Z 
pinch0’  but  having  higher  plasma  densities  such  as  those  of  the  gas-embedded  Z 

Q 

pinch.  In  addition  the  discharge  currents  in  the  plasma  focus  are  also  much 

higher,  which  reduces  finite  larmor  radius  (FLR)  effects.  Hence  understanding  the 

stability  of  the  plasma  focus  could  have  applications  to  other  types  of  Z  pinches 

especially  in  the  classical  MHD  limit  of  small  ratios  of  ion  gyroradius  (rcj)  and 

mean  freepath  ( to  plasma  radius  (ro;  i.e.  rcj/ro<0.1  and  X-Jr^OA). 

Besides  the  similarities  to  Z  pinch  approaches  to  controlled  fusion. 

understanding  the  plasma  focus  is  of  interest  in  its  own  right.  Until  1978  the 

9 

plasma  focus  held  the  record  for  the  most  fusion  neutrons  produced  per  shot.  It 
has  been  observed  in  large  plasma  focus  devices  that  the  majority  of  the  neutron 
production  occurs  during  the  stable  pinch  phase^;  hence  determining  the  physics 
behind  the  stability  of  the  plasma  focus  could  lead  to  a  viable  controlled  fusion 


concept.  Secondly,  soft  xrays  are  produced  during  the  pinch  phase  which  may  be 
used  in  soft  x-ray  lithography  and  microscopy11;  the  duration  of  the  soft  x— ray 
pulse  could  be  optimized  by  the  understanding  of  plasma  focus  stability.  Finally 
the  voltage  pulse  and  peak  current  interruption  generated  by  the  device  occur 
during  the  pinch  phase,1"  hence  investigating  the  stability  of  this  device  could  lead 
to  important  applications  in  pulsed  power. 

Experimentally  the  plasma  focus  discharge  is  observed  to  be  terminated  by 
a  Rayleigh— Taylor— like  m=0  (or  sausage)  instability.  However  the  period  of  the 
stable  phase  is  typically1  several  times  the  inverse  of  the  classical  growth  rate1'1 
for  this  instability.  In  this  paper  experimental  evidence  is  presented  for  an 
empirical  scaling  law  for  the  period  of  stability  or  more  directly  the  scaling  law 
for  h,  the  distance  along  the  plasma  column  between  the  site  of  the  m=0 
instability  and  the  end  of  the  column.  Then  the  theoretical  basis  for  such  scaling 
is  reviewed  and  these  scaling  laws  (both  empirical  and  theoretical)  are  applied  to 
predict  the  period  of  stability  in  other  Mather— type  1  plasma  focus  devices. 

In  the  next  section  the  phenomena  of  the  plasma  focus  discharge  related  to 
the  macroscopic  stability  (such  as  current  sheath  dynamics)  are  described.  In 
Section  III  the  experimental  apparatus  is  described  and  the  experimental  results 
governing  the  current  sheath  dynamics  and  the  scaling  of  h  and  rQ  with  plasma 
current  and  filling  pressure  are  presented.  In  Section  IV  the  theoretical  models 
are  formulated  to  predict  the  scaling  of  h  and  the  period  of  stability,  r.  with  r 
pinch  current  and  filling  pressure  and  compared  with  the  empirical  laws  presented 
in  Section  III.  Also  in  Section  IV  these  theoretical  and  empirical  laws  are  applied 
to  predict  r  in  other  Mather— type  plasma  focus  devices  and  compared  with 
experiment. 
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II.)  Relationship  of  h  to  stability 

After  the  current  sheath  of  a  Mather— type  plasma  reaches  the  end  of  the 
center  electrode  (Fig.  1)  the  current  sheath  collapses  radially  into  a  pinch  along 
the  device  axis.  On  reaching  the  axis  the  current  sheath  compresses,  radially 
expands  and  then  recompresses. 

IQ 

The  second  compression  is  associated  with  m=0  instabilities  which  occur 
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at  the  end  of  the  radial  expansion  phase.  ’  The  results  reported  here  are  in 
agreement  with  those  Findings.  The  predicted  period  of  stability  ^  for  a 
cylindrical  column  r  a  ^/Vg  is  about  5ns  (where  vg  is  the  velocity  of  sound) 
whereas  the  observed  stable  periods  are  typically  10  times  this.  For 

compressional  Z  pinches  the  enhancement  factor  is  at  least  20  which  is  usually 
explained  by  finite  Larmor  radius  (FLR)  effects.  But  due  to  the  higher  magnetic 
fields  present  in  the  plasma  focus  this  mechanism  should  be  ineffective. 

Alternatively,  the  pinch  column  is  observed  to  be  stable  as  long  as  the 
column  length  /  is  less  than  a  certain  number  a  of  column  diameters  d  where 

1  fv 

2<a<100  experimentally  and  most  typically  about  ten.  Since  the  current  sheath 
is  radially  concave  during  the  radial  implosion  phase,  the  pinch  first  occurs  in  the 
midplane  and  then  lengthens  as  the  trailing  portions  farther  and  farther  off  the 
mid— plane  reach  the  axis  in  the  so-called  zippering  effect  (Fig.  2).^  Thus  the 
apparent  speed  of  the  lengthening  of  the  column,  vc,  is  related  to  the  radial 
velocity  of  the  implosion,  vf  (which  is  roughly  a  constant  in  the  midplane,  Fig. 
3).  Therefore  one  might  expect  the  period  of  stability,  r,  to  scale  as: 
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where  f  is  a  constant  on  the  order  of  2  to  3  taken  from  experiment. 
Since  h  «  //2 


or 


oro 
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1 

so  if  v  a  v  a  v  then  one  would  expect  r  «  «'r  /v  ,  where  a'  =  0(7/2)  /-/f 
c  r  s  os 

for  8  =  1. 

To  illustrate  how  v£  and  vr  could  be  related  geometrically  consider  the 
following  simple  model.  Instead  of  a  cylinder,  consider  the  radially  collapsing 
current  sheath  as  a  concave  hyperbola  of  revolution  (Fig.  2)  moving  radially 
inward  with  a  uniform  v  .  As  the  first  compression  occurs  on  the  midplane  the 
surface  of  the  column  may  be  represented  by: 


72  ^  2 
“  kZ  +  ro 
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where  Z  is  the  distance  along  the  axis  from  the  midplane  and  r  is  the  radial 


position  of  the  current  sheath  at  height  Z  (Fig.  2). 

v  =k^v  where  v  =  Z.  Thus  one  would  expect  vc 
related,  where  f  =  in  equations  II  T.)  to  II  3.). 


For  Z,r>>rQ  one  has 

and  vr  to  be  directly 
This  zippering  effect 


occurs  in  all  Mather— type  plasma  focus  devices  and  has  been  reported  previously. 


Thus  this  elementary  model  predicts  h  should  be  directly  related  to  rQ  if  v 
and  vf  are  also  directly  related.  In  the  next  section  the  experimental  evidence  to 
support  this  model  will  be  examined. 


III.  Experimental  Procedures  and  Results 


Ilia.  Experimental  Apparatus  and  Diagnostic  Tool 

The  dense  plasma  focus  (DPF)  used  in  the  measurements  reported  here  is  a 
Mather-type1^  (Fig.  1).  The  center  electrode  is  hollow  and  has  an  outer  diameter 
of  5.08  cm.  the  inner  diameter  of  the  outer  electrode  is  10.2  cm  and  the  center 
electrode  extends  14  cm  beyond  a  lavite  insulator.  The  input  energy  on  the 
capacitor  bank  (40/iF)  was  varied  between  6  and  18  kj  in  these  measurements  and 
the  Filling  pressure  varied  between  1  and  8  torr  of  deuterium. 

The  velocities  of  the  plasma  edge  and  the  occurance  of  macroscopic 
instabilities  were  recorded  with  a  streak  camera  either  lined  up  with  a  vertical  slit 
to  observe  events  along  the  axis  of  the  pinch  (z  streak,  Fig.  4)  or  horizontally  to 
observe  the  radial  motion  of  the  plasma  column.  The  minimum  radius,  rQ  was 
also  taken  from  the  radial  streak  pictures,  whereas  h  was  taken  from  the  axial 
streak  pictures  (Fig.  5).  The  plasma  current  was  measured  with  magnetic-pickup 
loops  1mm  in  diameter  placed  in  a  quartz  tube  with  an  outer  diameter  of  3  mm. 
The  loops  were  calibrated  with  Helmhotz  coils  out  to  35MHz  and  placed  midway 
•between  the  inner  and  outer  electrodes  a  few  centimeters  below  the  top  of  the 
center  electrode.  Plasma  density  and  density  gradient  scale  lengths  of  the  pinched 
plasma  were  measured  for  some  discharges  by  use  of  a  Mach— Zehnder 
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interferometer  and  a  nitrogen  laser  (A=337nm). 
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Using  a  curved  crystal  spectrometer  ’  (Fig.  6),  the  electron  temperature 
of  our  DPF  was  determined  as  a  function  of  time.  After  a  line-free  region  of 
the  spectrum  was  located  (energy  regions  above  the  highest  line),  the 
bremsstrahlung  emission  was  monitored  with  three  PIN  diode  detectors.  This 
allows  the  determination  of  the  temperature  as  a  function  of  time  in  the  pinch. 
The  temperature  is  proportional  to  the  logarithmic  slope  of  the 


bremsstrahlung  spectrum  and  the  use  of  three  detectors  permits  the  verification  of 
the  effective  thermal  nature  (as  opposed  to  beam  target)  of  the  continuum. 

This  method  is  applicable  (with  our  high  efficiency  crystal  geometries)  even 
to  pure  hydrogen  or  deuterium  discharges.  There  is  no  need  to  use  a  seed  gas  to 
enhance  the  x— ray  emission,  unless  it  is  of  course  desired  to  investigate  the 
spectrum  of  such  seed  gas. 

For  hydrogen  and  deuterium  discharges  the  temperature  was  determined 
using  a  100  mm  x  20  mm  PET  crystal  of  eccentricity  0.806  placed  at  a  focal 
length  of  65  cm  from  the  Plasma  Focus  pinch.  The  energy  bands  selected  for 
this  application  were  1468  eV  (±5  eV)  and  1670  eV  (±6  eV).  The  spectrometer 
was  independently  calibrated  using  a  foil  filter  technique  (filtered  detectors  mounted 
at  the  end  of  a  4  m  collimating  tube)  and  a  second  crystal  spectrometer  (flat 
PET  crystals  at  1592  eV  and  2006  eV).  The  relative  efficiencies  of  the  PIN 
diode  detectors  were  also  determined  employing  the  long  collimating  tube  to  view 
the  same  region  of  the  plasma  pinch  during  a  series  of  calibration  discharges. 

The  temperature  obtained  using  this  technique  is  very  nearly  the  maximum 
temperature  obtained  at  the  most  dense  part  of  the  pinch,  as  it  developes  as  a 
function  of  time.  In  fact  the  linear  dependence  of  the  bremsstrahlung  radiation 
intensity  with  the  square  of  the  density  heavily  weighs  the  most  dense  regions  of 
the  pinch,  since  the  radiation  is  averaged  over  the  viewed  emitting  region. 

Ill  b.)  Experimental  Results 
Scaling  of  h  and  rQ 

The  experimental  scaling  and  magnitude  of  h  and  rQ  for  the  Illinois  DPF  as 
a  function  of  plasma  current  and  filling  pressure  are  shown  in  Figs.  7  to  10.  As 


can  be  seen,  the  power  law  least  squares  fits  to  the  data  reveal  that  h  depends 


on  the  plasma  current,  I,  and  filling  pressure  pQ  as: 


h .  ■  v 

where  x  =  — 0.02±0.12  and  y  =  0.26±0.13 


III  1) 


The  dependence  of  plasma  radius,  r  ,  with  these  quantities  as 


Ta  b 
r  oc  1  p^ 
o  Fo 

where  a  =  -0.18±0.11  and  b  =  0.2±0.03 


III  2) 


Thus  within  the  uncertainties  h  and  rQ  have  the  same  dependence  on  these 

quantities  and  hence  are  proportional  to  each  other  within  the  experimental  error. 

There  is  a  shift  between  the  current  scaling  value  for  r  with  r  = 

1.7  ±  0.2  mm,  and  the  pressure  scaling  value  with  rQ  =  1.26  ±  0.8mm.  Although 

within  experimental  uncertainties,  the  discrepancy  may  be  due  to  slight  differences 

in  the  camera  alignment  since  these  data  sets  were  taken  on  two  separate 

occasions.  A  typical  value  of  rQ  obtained  from  the  Mach-Zehnder  interferometer 
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(half— width  half— maximum  of  the  density  profile)  is  1.5±0.1mm  (Fig.  11),  which 
is  within  the  uncertainties  of  both  sets  of  dates  obtained  from  the  streak  camera. 


Relationship  Between  v  and  v 

r  c 

The  streak  camera  was  lined  up  horizontally  to  view  the  region  0.27±0.05 
cm  above  the  top  of  the  center  electrode.  The  current  flowing  in  the  plasma  was 
monitored  by  a  magnetic— pickup  loop  whose  center  was  2.4  cm  below  the  top  of 
the  center  electrode  and  3.4  cm  from  the  device  axis.  Typical  results  are 
presented  in  Fig.  3  showing  that  the  radial  collapse  velocity  achieves  a  roughly 


10 


constant  value.  The  results  of  the  scaling  of  the  radial  velocity  at  a  given  radial 

position  as  a  function  of  plasma  current  and  filling  pressure  are  shown  in  Figs.  12 

and  13.  For  the  scaling  of  vf  with  current,  the  prediction  of  the  "M"  theory 
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adapted  to  the  Illinois  DPF  is  also  shown. 

The  results  versus  current  do  not  go  through  zero  as  predicted  by  theory, 

23 

according  to  which  vf  should  scale  as: 
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UtA  (2"omino> 
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where  v^  is  the  Alfven  speed  at  the  outer  edge  of  the  center  electrode  (r=r^),  I 
is  the  plasma  current,  nQ  is  the  fill  density  and  mj  is  the  ion  mass.  The 
deviation  of  the  data  from  the  theory  may  be  due  to  a  velocity  dependent  drag  of 
the  current  sheath  on  the  anode  but  this  is  only  speculation  on  our  part. 

vc  was  measured  from  the  axial  streak  pictures  (Fig.  5)  and  is  plotted  as  a 
function  of  I  and  pQ  in  Figs.  14  and  15.  To  see  if  vr  and  v£  are  related  we 
compare  the  least  square  fits  to  the  data  presented  in  Figs.  12  to  15. 

For  the  dependence  of  vr  and  vc  with  plasma  current  the  two  data  sets 
were  fitted  to  a  power  law  relationship,  i.e. 

v.  =  B  IR  III  4.) 

and  vc  =  BpI^  III  5.) 

where  R  =  0. V 1  ±0.07,  and  F  =  0.78±0.07 

Thus  within  the  uncertainties  (less  than  10%)  these  two  velocities  have  the  same 
dependence  on  the  plasma  current  as  would  be  expected  in  the  zippering  effect. 
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For  the  variation  in  pressure  vr  and  vc  were  fit  to  a  power  law: 


vr  «  Po 


and  vc  «  p0 


where 


III  6.) 

Ill  7.) 


r  =  — 0.41±0.04,  and  c  =  — 0.46±0.04 


Thus  within  the  10%  experimental  uncertainties  vf  and  vc  have  the  same 
dependence  on  filling  pressure. 


Plasma  Electron  Temperature 

The  electron  temperature  is  useful  in  calculating  the  speed  of  sound  in  the 

plasma  for  comparison  with  theory.  The  soft  x— ray  measurements  involve  an 

average  over  a  few  millimeters  along  the  pinch  column.  For  pure  hydrogen  or 
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deuterium  filling  gases  (Fig.  16),  the  electron  temperatures  typically  rises  to 
250eV  at  the  peak  of  first  compression,  cools  to  about  200eV  at  the  end  of  the 
radial  expansion  (about  16ns  later)  and  rises  back  up  to  over  300eV  during  the 
peak  of  second  compression  (about  27ns  after  first  compression). 

The  results  presented  in  this  section  will  be  discussed  in  the  following 
section. 


IV)  Discussion  of  Results 

In  this  section  the  experimental  results  are  summarized  and  then  compared 
with  theory.  Since  the  stability  period  appears  to  be  enhanced,  two  simple  models 
are  discussed  which  could  account  for  such  stability  and  still  scale  as 


I<i  k'i'il|'tl'.V«<<  i<l  ki  k«# 


rQ/vg.  Finally  the  r  a  a'  JxQ  scaling  is  compared  to  data  from  other  devices  to 
check  its  validity. 

IVa)  Summary  of  Experimental  Results 

From  the  discussion  of  section  II  it  is  seen  that  is  h  «  rQ  and  that  if 
vc  oc  vf,  the  scaling  would  be  consistent  with  the  classical  result  that  t  s  ‘  a 
rQ/vg.  h  (the  distance  between  the  point  the  m=0  stability  occurs  and  the  end 
of  the  pinched  column)  is  a  convenient  parameter  to  measure  from  axial  streak- 
pictures.  The  experimental  test  of  these  relationships  was  to  vary  the  plasma 
current,  I  and  filling  pressure  (deuterium),  pQ,  over  as  wide  a  parametric  range 
possible  to  see  how  each  of  these  quantities  varied  with  these  parameters. 

For  variations  of  the  pair  h  and  rQ  with  pQ,  and  for  variations  of  the  pair 
vf  an  vc  with  both  pQ  and  I  the  same  scaling  is  found  (see  egns.  Ill  1  to  III  7). 

The  vf  and  vc  scaling  with  pQ  and  I  give  support  to  the  zippering  model 
of  the  relation  between  these  two  quantities,  h  scales  the  same  way  as  rQ  with  I 
within  the  uncertainties  which  are  large  with  respect  to  the  mean  value  (egn.  Ill 
1  and  2)  but  which  are  at  most  less  than  twice  those  of  the  vr,  vc  data  in  the 
absolute  sense.  Since  the  mean  values  of  the  exponents  are  close  to  zero  this 
merely  implies  that  the  dependence  of  h  and  rQ  on  I  is  weak.  While  this  does 
not  prove  that  h  and  rQ  are  related,  the  demonstrated  correlation  between  vc  and 
vf  supports  the  classical  scaling  model: 


IV  1) 


and  the  weak  dependence  of  h  and  rQ  with  I  is  at  least  consistent  with  this. 
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Also  the  small  but  significant  trend  in  the  scaling  of  h  and  rQ  with  pQ  (i.e.  h  * 
pQ  0.26±0.13  and  rQ  «  pQ  0.2±0.03)  indicates  that  h  and  rQ  are  directly  related. 

a,  the  ratio  h/rQ,  is  8.6±2.4  experimentally  and  related  to  a'  =  rQ  by 
(section  II) 

a'  =  o(7/2)1/2/f.  IV  2.) 

where  f  =  v  /v  f  is  2.6±0.8  for  p  variations  and  1.4±0.7  for  I  variations 
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making  the  mean  value  f  a  1.9±0.9  .  Thus  a'  *  4  ±  2  from  the  experimental 
results  of  h/rQ  and  vc/vr  and  this  model. 

This  will  be  compared  with  a'  taken  from  other  experimentally  determined 
quantities  and  theoretical  models  for  a '  in  the  next  subsection. 

IV  b)  Comparison  of  Experiment  with  Theoretical  Models  of  the  Sausage 

Instability  in  the  Plasma  Focus 

In  this  section  the  classical  Kruskal-Schwarzschild3  theory  for  the  sausage 
mode  is  applied  to  the  pinch  phase  of  the  plasma  focus  and  used  to  predict  the 
growth  rate  and  wavelength  of  the  instability.  Then  three  additional  models  for 
enhanced  stability  based  on  1)  the  shape  of  the  pinch,  2)  axial  flow  of  plasma, 
and  3)  self-generated  magnetic  fields  are  also  presented. 

The  Kruskal-Schwarzschild  Theory 

The  classical  mhd  theory  for  the  m=0  mode  in  a  plasma  supported  by  a 
magnetic  field  against  gravity;  this  is  the  well-known  Kruskal-Schwarzschild 
instability.  If  the  gravitational  acceleration  vector,  g,  points  away  from  the  fluid 
the  plasma  is  unstable  and  the  growth  rate  of  the  instability,  f,  is  related  to  the 
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wave  number  k  by: 


r  =  (kg)1/'2  IV. 3) 

for  the  case  k  •  B  =0.  The  gravitational  acceleration  vector  pointing  out  of 
the  plasma  is  equivalent  to  the  plasma  being  accelerated  by  the  magnetic  field  or 

a  convex  curvature  of  the  magnetic  field  with  respect  to  the  plasma;  the  latter 
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causes  the  centripetal  acceleration  vector  of  the  fluid  to  point  into  the  plasma.-' 
Thus  during  first  compression,  the  plasma  is  being  decelerated  by  its  own 
compression,  the  acceleration  vector  points  out  of  the  plasma,  and  hence  is 
stable  .  However  toward  the  end  of  the  radial  expansion  phase  the  plasma  is 
being  decelerated  by  the  magnetic  field  and  hence  should  be  unstable. 

Although  equation  IV.2  predicts  an  infinite  growth  rate  as  the  wavelength  of 
the  perturbation  approaches  zero,  two  phenomena:  finite  ion  gyroradius,  p  and 
finite  density  scale  length,  77  \  can  give  a  lower  bound  to  the  wavelength.  When 
one  considers  the  effect  of  field  line  curvature  for  a  cylinder  g  can  be  given  by:*1 

g  =  «T/mjr0  IV.4) 

(here  T,  mj  and  rQ  are  the  ion  temperature,  ions  mass  and  the  plasma  radius 

27 _ 28 

respectively),  one  finds: 


From  our  experiment  during  first  compression,  where  pQ  =  3  torr  (D9),  I  = 
550kA,  T~  330eV2^,  rQ  ~  1.5mm  and  t)  ~  (1.2±0.2mm)-1  (Fig.  11)  r  ~  1/6  - 
ll±2ns.  This  compares  with  the  experiments  period  between  first  and  second 
compression  of  27±3ns  so  a'  for  our  experiment  from  the  classical  model  is  2.5±0.5 
which  is  consistent  with  the  data  from  the  hA0>vc>vr  scaling  results. 

IV  c)  Effects  of  Curvature 

During  the  radial  collapse  phase  the  plasma  experiences  an  acceleration  just 
prior  to  first  compression  (Fig.  3)  which  could  cause  m  =  0  modes.  While 
peturbations  in  the  current  sheath  are  commonly  observed  during  this  period 
(indeed  there  is  evidence  of  such  in  the  upper  right  portion  of  the  interferogram 
in  Fig.  11),  they  don't  cause  a  column  disruption  as  in  second  compression.  Since 
the  column  is  concave  during  radial  collapse  it  is  possible  that  the  growth  of 
these  perturbations  is  surpressed  by  curvature  effects  considered  by  Comisar 
Later,  during  the  radial  expansion  phase,  the  plasma  column  is  cylindrical  near  the 
midplane,  so  the  effect  of  curvature  would  be  reduced.  Therefore  some  other 

mechanism  must  provide  stability  during  this  phase. 

IV  d.)  Stability  induced  bv  Axial  Flow 

During  first  compression  the  plasma  is  stabilized  by  its  own  deceleration  but 
during  the  radial  expansion  this  mechanism  is  no  longer  effective.  However  there 
is  considerable  evidence  of  axial  fluid  flow  during  this  phase  which  could  provide 
stabilizing  effects. 

Besides  the  radial  motion  of  the  plasma  during  these  phases,  strong  axial 
flow  is  also  implicit  due  to  the  lack  of  confinement  at  the  ends  of  the  pinch. 

91 

Downing  and  Eisner  using  laser  scattering  during  the  pinch  of  a  15kJ 


Mather— type  plasma  focus  observed  a  roughly  constant  density  during  the  last 
stages  of  the  radial  collapse  through  first  compression.  Since  the  radial  velocity 

*■*  o 

was  ~3*10<  cm/sec  during  this  period,  an  axial  velocity  of  2.4.10  cm/'sec  was 

implied  to  conserve  particles.  They  also  analyzed  other  laser  scattering 

observations  of  the  plasma  focus  and  concluded  high  axial  velocities  were 

probably  present  in  these  devices  as  well.  Also  there  is  evidence  of  high  velocity 

(~10°  cm/sec)  axial  flow  from  x— ray  framing  camera  pictures.  The  presence  of 

a  high  axial  velocities  during  the  first  compression  and  radial  expansion  phases  is 

35—37 

supported  by  computer  simulations  of  these  phases  as  well. 

In  our  device,  evidence  of  high  speed  axial  flow  can  be  seen  in  the  rapid 
decrease  in  particle  number  per  unit  length,  N,  during  the  three  phases  of  the 
pinch  as  measured  using  the  Mach-Zehnder  interferometer.  The  interferograms 
reveal  that  for  the  standard  case:  pQ  =  3  torr  deuterium  and  I  =  550k A,  N  falls 

10  7  7/ 

from  about  10  /cm  at  first  compression  to  ~  6.10  'cm  during  the  radial 
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expansion  and  down  to  10  /cm  during  second  compression.  Thus  the  rate  of 
decay  of  N  is  roughly  constant  of  over  this  period. 

Taking  a  simple  axial  outflow  model: 

/dN  =  2Nv7  IV  6.) 

d  t  Z 

where  /  v  are  the  length  of  the  pinch  area  and  the  axial  flow  velocity 

respectively.  Since  /  =  2h  and  taking  an  average  h  ~  1cm  during  radial 

7 

expansion,  one  finds  v z  ~  6.10  cm/sec  for  the  entire  pinch  period. 

OO 

The  stability  critrion  may  be  expressed  by  v^  >  v^;  perturbations  are 
swept  at  out  of  the  column  before  they  can  grow.  However  as  the  mass  flows 


out  the  ends,  increases  until  the  instability  is  triggered.  For  I  ~  550kA  and 

7 

pQ  ~  3  torr,  analysis  of  the  interferograms  reveals  that  ~  4.6  •  10  cm/sec 
at  the  time  of  second  compression.  The  speed  of  sound  on  the  other  hand,  since 
200eV,  is  ~1.4  •  10 1  cm/sec  through  out  the  pinch  phase,  thus  the  flow  is 
supersonic  with  a  Mach  number  of  ~  4. 

In  this  model  the  stable  period  scales  as  a  loss  time  for  mass 

T{  =  h/vZ 

Theoretically  can  be  related  to  vr  as  follows.  Because  of  particle 
conservation  during  the  radial  collapse,  the  axial  flow  velocity  vz  for  a  cylinder  of 
length  /  and  radius  r  would  be  ~  /  vr/r  where  vr  is  the  radial  velocity.  This 
flow  velocity  appears  to  persist  well  into  the  radial  expansion  phase.  This  may 
be  seen  theoretically  as  follows:  If  one  assumes  potential  flow,  then  Bernoulli’s 
equation  holds  and  one  has: 


p/p  +  1/2  v  £  =  constant 


IV  7.) 


The  first  term  in  IV  7.)  is  just  v:/ 7,  where  7  is  the  ratio  of  specific  heats. 

b 

Since,  as  shown  above,  v  <<v  ),  v7  would  be  approximately  constant  during  the 
entire  pinch  as  consistent  with  the  above  discussion. 

From  these  considerations  appears  that  v^  <  vz  during  radial  expansion, 
but  this  reverses  at  2nd  compression  due  to  mass  loss.  The  scaling,  since  v  « 
vr  would  be  essentially  classical: 


t  ~  h/v 
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with  about  the  enhancement  factor,  o'  =  o/  /  (7/2)^“  being  somewhat  smaller 

r 

(~2.6)  since  v7  >  v  .  Thus  this  model  also  appears  to  give  results  consistent 

la  C 

with  the  data. 


IV  e)  Enhanced  Stability  from  a  Self  Generated  Magnetic  Field 

A  self  generated  axial  magnetic  field  could  provide  column  stability.  The 
addition  of  an  external  axial  magnetic  field  has  been  shown  to  provide  stability  to 

QQ 

the  plasma  focus;  the  elementary  theory  of  the  Kruskal— Schwarzschild 

instability25,4®  also  supports  this.  Filamentary  structures, observed  in  the  discharge 
from  the  beginning  ’  “  and  found  in  the  pinch  could  be  force— free 
configurations44  hence  providing  an  axial  magnetic  field.  Evidence  for  such  a  field 
has  been  reported  from  measurements  using  magnetic  pick-up  loops41,45  and 


Faraday  rotation. 


However  in  ref.  41  the  filamentary  fields  were  found  to 


alternate  between  inwardly  and  outwardly  pointing  radial  fields  during  the  rundown 

phase  whereas  in  ref.  45  they  were  found  to  point  only  in  one  direction  radially. 

Also  in  ref.  42  the  filamentary  structures  were  reported  to  disappear  later  in  the 

discharge  whereas  in  ref.  41  they  were  reported  to  persist  into  all  phases.  Thus 

the  nature  of  filamentary  structures  is  somewhat  ambiguous. 

A  possible  model  to  take  into  account  the  magnetic  shear  introduced  by  the 

axial  magnetic  field  of  such  filaments  could  be  as  follows.  The  shear  limits  the 

25 

wave  number  of  the  maximun  gr  /th  rate  ,  k^^  t0 


MAX 


kT  popo 
mi  ro  2  B  2 


IV  9. 
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where  Bz  is  the  axial  magnetic  field.  Assuming  (1  =  1,  the  wavelength  of  the 
maximum  growth  rate,  'r  now  S’ven  by; 


/-  2h>AMAX/2 


IV  11.) 


or 


h> 


4,ro 


IV  12.) 


For  comparison  with  our  experiment,  o=h/rQ  ~  8.6,  Bz/ ~  0.6.  The 

magnitudes  of  the  Br  and  Bz  and  B^  magnetic  fields  reported  in  ref.  41.45  and 

46  would  make  Bz/B^  of  this  order  or  greater  so  this  effect  cannot  be  ruled  out. 

However,  the  observed  wavelength  of  the  instability  is  at  least  an  order  of 

magnitude  smaller  than  X  =  4h  as  predicted  by  the  model,  possibly  due  to  its 

simplicity.  This  model  is  only  presented  here  as  an  illustration  of  one  effect 

44 

filamentary  structures  could  have  on  stability.  A  more  sophisticated  treatment 
is  beyond  the  scope  of  this  work. 

Computer  simulations  of  the  plasma  focus  dynamics  using  a  two— fluid  model 
reveal^"  that  current  filamentation  should  occur  in  the  pinch  phase  which  could 
explain  the  structures  observed  by  Sadowski  et  al.^  and  provide  a  stabilizing  Bz 
as  well.  The  results  of  these  simulations  also  show  that  the  period  of  stability 


should  scale  as  the  plasma  current,  I,  to  the  three— halves  power.  Because  of 
the  complicated  nature  of  these  calculations,  it  is  not  obvious  that  there  is  a 
single  mechanism  involved  in  the  enhanced  stability. 


Mather— tvne  Plasma  Focus  Devices. 

So  far  evidence  has  been  presented  to  indicate  that  h  =  a  rQ  and 

t  =  <*'r0/vr  in  our  Mather-type  plasma  focus  but  is  this  result  unique  to  this 

device?  Since  the  scaling  is  classical,  enhanced  by  a  factor  of  a',  it  seems 

reasonable  that  it  should  apply  to  other  DPFs.  We  have  gathered  data  or  pinch 

current,  fill  pressure  (i.e.  equivalent  in  mass  density  of  that  of  pure  deuterium), 

pinch  radius,  and  r  from  the  literature  and  scaled  v  «  j0.745^  q  435  where  the 

r  Po 

exponents  are  the  mean  values  from  our  v  vf  data  to  be  consistent  with  the 

zippering  effect. 

The  results  are  shown  in  Table  IV  1.)  and  Fig.  17.  The  scaling  relations 

for  the  stability  period  is  given  by 


„  ro  r1  .10.745  po  lo.435 

7\  =  T  . 

5  4  rn  ll  Lp_. 


IV  13.) 


where  the  quantities  subscripted  by  c'  are  those  for  our  device  (U.  of  I.)  and  are 

listed  in  Table  VI  1.)  along  with  those  parameters  from  the  other  devices.  The 
data  follow  the  linear-least  square  fit  fairly  well  with  less  than  a  20%  uncertainty 
in  the  slope  (1.3  ±  0.2).  All  those  listed  but  the  Frascati  device  have  hollow 

cathodes  which  could  effect  vf  if  drag  on  the  center  electrode  is  an  important 
effect.  That  is,  with  a  hollow  cathode,  the  drag  force  is  farther  away  from  the 

constricting  region  during  collapse  than  with  a  solid  anode  and  hence  the  effects  of 
drag  would  be  reduced.  This  could  explain  the  devation  of  the  Frascati  point  (F) 
from  the  rx  =  rg  line  but  not  that  of  the  Poseidon  point  (P).  In  any  case  the 

agreement  with  a  linear  relationship  between  r  and  r  is  evidence  that  the 
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scaling  relation  IV  13.)  is  on  the  right  track.  The  values  of  hy  were  scaled 
simply  from: 


h 


s 


IV  14.) 


however  fewer  values  of  h^  can  be  found  in  the  literature  (Table  IV  1). 
However,  those  that  are  found  are  in  fair  agreement  with  the  predictions  of  IV 
14.  Thus  this  enhanced  classical  scaling  appears  to  hold  for  several  other 
Mather— type  plasma  focus  devices  as  well. 

Although  the  scaling  equations  IV  13.)  and  IV  14.  are  useful  to  test 
consistency  between  devices,  they  are  not  very  useful  in  design  since  one  has  to 
know  rQ  before  one  can  predict  r.  If,  however,  one  plots  rQ  versus  rA,  the  anode 
radius,  for  the  same  devices  of  Table  IV  1.)  one  finds  (Fig.  18): 

rQ  =  (6±  1.6)  •  1(T2  rAL3±0-2  IV  15.) 


both  rQ,rA  in  cm. 

If  the  radial  collapse  phase  were  adiabatic  and  one-dimensional,  one  would 
2 

predict:  (since  V  «  r  ) 


Since  it  has  been  noted  that  the  pinch  temperatures  TQ  do  not  change 

55 

significantly  over  a  wide  range  of  devices  ,  the  near  linearity  between  rQ  and  r  ^ 
is  not  too  surprising.  ' 

Thus  it  appears  one  could  use  the  scaling  relation  IV  15.)  along  with  IV 
13.)  and  IV  14.)  to  predict  r  and  h  for  a  given  design.  However  due  to  the 
scatter  in  the  data,  (the  uncertainties  in  these  predictions  approach  50  percent) 
only  rough  estimates  can  be  obtained  and  the  relations  strictly  apply  only  to 
devices  with  a  hollow  anode. 

Finally  the  empiracle  scaling  of  r  with  I  from  egn.  IV  13.)  would  appear  to 
be  at  odds  with  the  four— fluid  model  developed  at  the  University  of  Stuttgart 
which  predicts47  r  «  I3^.  However,  if  rQ  is  plotted  verses  I  for  the  various 
devices  listed  in  table  IV  2.)  one  finds  that: 

r  *  I  L44±0-2  IV  17.) 


from  a  least-square  power  law  fit  so  the  overall  empiracle  scaling  is: 


T 


a  r  t0.7±0.3  0.435±.04 

°  «  !  pQ 


IV  18.) 


While  this  is  still  not  within  the  experimental  uncertainties  of  the  computational 
results47,  the  disagreement  is  not  nearly  as  severe  as  IV  13.)  would  indicate. 
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Table  IV  1 

Scaling  of  h  and/or  r  to  other  Mather-type 
Plasma  Experiments  Focus  Devices:  Comparison 
with  Experiments 


Exp 

Exp 

Empirical 

Scaling 

Location  of 

rA 

I 

Po 

ro 

T 

h 

Device 

(cm) 

(kA) 

(torr) 

(cm) 

(ns) 

(cm) 

r(ns)  h(cm 

U.  of  Illinois 

2.5 

550 

3.0 

0.15 

27 

1.3 

27  1.3 

14  4fi 

Frascati  ’ 

8.0 

1300 

3.0 

1.0 

130 

— 

95  8.7 

Stuttgart 

Poseidon10’45’48 

6.7 

1200 

2.3 

0.78 

100 

— 

70  6.S 

Nessi'y 

3.3 

290 

1.5 

0.15 

40-60 

2 

32  1.3 

Swierk50^ 

5.0 

1000 

II 

0.38 

75 

2-2.5 

71  3.3 

Darmstadt52’53 

0.8 

200 

2.6* 

0.04 

16 

0.25 

14.4  0.35 

za 

Stevens 

1.7 

400-500 

8 

0.2 

60±20 

1.0-1. 5 

64+6  1.7 

Equivalent  pressure  in  deuterium  of  4  torr  Hg  +  4%  Ar  in  terms  of  mass 
density 
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Earlier  data  reported  from  Swierk  gave  a  somewhat  greater  plasma  radius 
as  measured  using  a  streak  camera.  The  data  reported  in  ref.  50  are  taken 
from  Schlieren  pictures  and  laser  interferograms  and  hence  are  more  accurate. 


VI  Conclusions 


The  experimentally  observed  constancy  of  the  ratio  of  the  distance  between 
first  and  second  compression  along  the  axis  of  the  pinch,  h,  and  the  minimum 
plasma  radius,  r  ,  is  evidence  that  the  stability  of  the  Mather— type  plasma  focus 

is  behaving  in  a  manner  as  predicted  by  ideal  MHD  theory.  The  scaling  of  these 

parameters  with  plasma  current  and  filling  pressure  are  reported  here  for  the  first 
time.  Although  the  ratio  a  =  h/rQ  is  about  a  factor  of  three  above  that 
predicted  for  a  straight  cylinder,  two  simple  models,  both  consistent  with  ideal 

MHD  theory,  can  be  used  to  explain  this  enhancement.  The  models  include  either 
the  effect  of  axial  flow,  or  the  effect  of  a  self-generated  axial  magnetic  field. 

However  neither  of  these  models  is  necessarily  unique. 

The  empirical  scaling  laws  developed  from  the  experimental  results  found  in 
this  experiment  were  used  to  predict  both  h  and  the  period  of  stability  r  in  other 
plasma  focus  devices.  The  agreement  between  the  scaled  parameters  and  those 
reported  from  the  other  experiments  (table  IV  1)  indicates  that  the  empirical 
scaling  relations  for  h  and  r  foumd— tn  this  study  could  be  used  to  predict  the 
performance  of  other  devices.  Thus  if  enhanced  stable  lifetime  of  the  pinch  is 
desired  at  higher  plasma  currents,  one  could  design  a  device  with  an  anode  of 
larger  radius  (egns.  IV  13—15). 
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Fig.  1  Schematic  of  a  Mather— type  plasma  focus  showing  the  successive 

positions  of  the  current  sheet  and  the  pinch  plasma  at  the  end  of  the 
center  electrode. 

Fig.  2  The  relationship  between  the  radial  velocity  and  the  apparent  vertical 

speed  of  first  compression.  If  the  shape  of  the  current  sheath  and  the 
radial  velocity  prior  the  arrival  on  the  axis  are  held  constant  in  space 
and  time  first  compression  would  appear  to  move  along  the  axis  as  in 
b).  The  relationship  between  radius  and  height  of  the  current  sheath 
is  taken  to  be  hyperbolic. 
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Fig.  3  A  comparison  between  the  predictions  of  "M"  theory  ’  and 

experimental  observations  using  radial  streak  pictures  and  magnetic 
pickup  loops.  The  width  of  the  vertical  bars  on  the  data  are  the 
experimental  ranges  of  uncertainty. 

Fig.  4  Alignment  of  the  slit  for  the  streak  camera  for  the  z  streak  pictures. 

Fig.  5  Examples  of  z  streak  camera  pictures  (top  and  middle)  and  the 

corresponding  dl/dt  waveform,  h  the  distance  between  first  and 
second  compression  when  second  compression  occurs  is  defined  in  t.he 
middle  figure. 

Fig.  6  Elliptical  curved  crystal  spectrometer  used  to  determine  the  electron 

temperature.  The  plasma  column  is  at  one  focus  of  the  ellipse  and  the 
cross  over  point  at  the  other.  Pin  diodes  on  the  detection  circle  are 
used  to  detect  different  spectral  bands  of  the  soft  x— ray  emisson. 

Fig.  7  The  dependence  of  h  (Fig.  5)  on  plasma  current.  The  deuterium 

filling  pressure  was  3.0  torr. 

Fig.  8  The  dependence  of  h  (Fig.  5)  on  deuterium  filling  pressure.  The 

initial  energy  on  the  capactior  bank  was  12kJ. 

Fig.  9  The  dependence  of  the  minimum  plasma  radius,  r  ,  on  plasma 

current.  The  deuterium  filling  pressure  was  3.0  torr. 

Fig.  10  The  dependence  of  the  minimum  plasma  radius,  r  ,  on  plasma 

current.  The  initial  energy  on  the  capacitor  bank  was  12kJ. 

Fig.  11  Top:  Typical  interferogram  of  the  first  compression  phase  obtained 
from  the  nitrogen  laser  interferometer  (A=337  nm,  Ar<5ns).  Bottom: 
Abel  inversion  of  interferogram  yields  this  electron  density 
20 

distribution. 


Fig.  12  The  radial  velocity  during  the  collapse  phase  at  a  radius  of  1.1  cm  and 
a  filling  pressure  of  3  torr  D9  as  a  function  of  plasma  current.  The 

dashed  line  and  solid  line  are  the  least— squares  power-law  fit  to  the 

9Q  94 

data  and  the  predictions  of  "M"  theory"  respectively. 

Fig.  13  The  radial  velocity  during  the  collapse  phase  at  a  radius  of  0.8cm  and 
a  bank  energy  of  12kJ  as  a  function  of  deuterium  filling  pressure. 

Fig.  14  Apparent  axial  speed  of  first  compression  about  1  cm  above  the 
midplane  as  a  function  of  plasma  current.  The  deuterium  filling 
pressure  was  3.0  torr.  The  line  is  a  least  squares  power  fit  to  the 
data. 

Fig.  15  Apparent  axial  speed  of  first  compression  about  1  cm  above  the 
midplane  as  a  function  of  deuterium  filling  pressure.  The  initial 
energy  on  the  capacitor  bank  was  12kJ. 

Fig.  16  Typical  time  dependence  of  the  electron  temperature  of  the  pinched 
plasma  obtained  with  the  curve  crystal  spectrometer  for  P0  =  3  torr 

(deuterium)  and  I  =  550k A.  Comparison  with  simultaneous  laser 

interferograms"  indicates  the  first  peak  corresponds  to  first 
compression,  the  dip  to  radial  expansion  and  the 
second  peak  to  second  compression. 

Fig.  17  A  comparison  of  the  experimental  period  of  stability,  r  ,  with  that 

predicted  from  egn.  IV  16.).  The  solid  line  is  the  least  squares  linear 
fit  to  the  data  Key:  F  =  Frascati,  P  =  Stuttgart  Poseidon,  N  = 
Stuttgart  Nessi,  D  =  Darmstadt,  St  =  Stevens  Inst.,  Sw  =  Swierk, 
and  I  =  Univ.  of  Illinois  (Table  IV— 1). 

A  plot  of  the  final  pinch  radius,  rQ  versus  the  radius  of  the  anode,  r^, 

for  the  devices  listed  in  Table  IV-1.  The  line  is  the  least  squares 
power  law  fit  to  the  data. 
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I.  INTRODUCTION 

The  determination  of  electron  density  in  a  plasma 
is  important  in  the  formulation  of  plasma  models.  Some  of 
the  more  common  parameters  that  the  electron  density,  ne, 
helps  to  reveal  are  electron  plasma  frequency  (fp  a  ne^), 
electron  collision  rates  CR  a  ne),  plasma  skin  depth 
(£s  a  ne-^),  and  Debye  length  (\D  a  ne-^). 

The  time  resolved  density  profile  is  important  in 
displaying  plasma  boundaries  and,  resolution  permitting, 
finer  structures  of  the  plasma  such  as  current  sheaths.  In  a 
dynamic  plasma  such  as  that  generated  in  a  plasma  focus, 
knowledge  of  the  density  profile  is  fundamental  to  the 
understanding  of  plasma  behavior  and  permits  comparison  with 
various  theoretical  models. 

It  is  the  purpose  of  this  work  to  measure  the 
electron  density  profile  of  the  plasma  produced  by  a  dense 
plasma  focus  in  order  to  study  the  dynamics  of  the  pinch 
phase;  a  nitrogen  laser  interferometer  is  used. 

A.  The  Dense  Plasma  Focus 

The  dense  plasma  focus  CDPF)  is  a  machine  that 
transforms  a  cold,  low  pressure  gas  into  a  hot,  high  density 
plasma.  Conversion  is  achieved  by  the  magnetic  pinch  effect. 
Plasma  densities  on  the  order  of  lO^cm-3  are  produced  along 
with  high  neutron  yields  and  intense  radiation  across  a  wide 
region  of  the  spectrum- 
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Because  significant  nurr;v  v.  of  neutrons  are  pro¬ 
duced,  the  DPF  was  originally  stud;-  -s  a  potential  fusion 
reactor.  However,  the  dense  plasma  c  not  last  long  enough 
in  present  day  devices  to  obtain  .  "creak  even"  reactor-*-. 
The  DPF  is  still  being  researched  as  a  high  energy  opening 
switch^,  soft  x-ray  source^,  laser  pumping  device4,  neutron 
source  and,  aside  from  many  other  possible  applications,  it 
is  being  studied  for  a  better  understanding  of  hot,  dense 


plasmas. 


The  device  consists  of  a  oair  of  coaxial  electrodes 


-  open  at  one  end  and  insulatec  a;  the  other  -  in  a  gas 
filled  chamber  (Figure  1).  Up  to  100  kJ  of  electrical  energy 
stored  in  a  capacitor  bank  is  discharged  into  the  electrodes 
through  spark  gap  switches.  The  pinch  process  starts  as  an 
electrical  breakdown  of  the  gas  along  the  base  insulator 
forming  a  current  sheath  between  the  electrodes.  The  self- 
induced  magnetic  field  lifts  the  current  sheath  from  the 
insulator  and  accelerates  it  down  the  electrodes.  When  the 
open  end  of  the  electrodes  is  reached,  the  current  sheath 
collapses  radially,  forming  a  high  density  plasma  column. 

For  the  University  of  Illinois  DPF,  the  center 
electrode  has  a  diameter  d^=5  cm,  the  outer  electrode  has  an 
inner  diameter  dQ=10  cm  and  the  length  is  25  cm.  The  energy 
is  initially  stored  in  four  10  mF  ■_  apacitor  banks  placed 
symmetrical ly  around  the  electrodes-  ~ne  initial  voltage  on 
the  capacitors,  VQ,  can  be  varied  between  20  kV  and  30  kV  in 
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order  to  vary  the  plasma  current.  The  gas  in  the  fill 
chamber  is  usually  deuterium  or  hydrogen  at  3  torr.  The 
pressure,  PD,  can  also  be  varied  from  1.5  torr  to  5  torr. 

B.  Overview  of  Interferometry  as  a  Plasma  Diagnostic 

All  interferometers  rely  on  interference  to  produce 
fringes  on  a  screen  (Figure  2).  A  single  electromagnetic 
wave  is  split  in  two  and  recombined  with  a  slight  angle,  0, 
between  the  wave  vectors,  and  k2-  The  phase  fronts 
interfere  constructively  and  destructively  producing  light 
and  dark  areas  on  the  screen  known  as  fringes.  For  a 
rigorous  mathematical  treatment  of  fringe  production  see 
Reference  5. 

The  two  waves  must  originate  from  a  single  source 
in  order  to  preserve  the  phase  information.  It  is  easily 
seen  from  Figure  2  that  the  fringes  will  move  across  the 
screen  proportionally  with  a  change  in  phase  between  the  two 
waves.  The  phase  shifts  can  be  random,  caused  by  the  natural 
line  broadening  of  the  source,  or  determined  by  a  change  in 
optical  path  length  of  one  wave  with  respect  to  the  other. 
The  second  case  is  of  interest  here.  If  a  material  of  length 
L  and  index  of  refraction  pCy)  is  inserted  in  one  arm  of  the 
interferometer,  the  fringes  will  change  position  on  the 


screen  according  to 


where  6  is  the  fringe  shift  in  wavelength  units,  X  is  the 
source  wavelength  in  a  vacuum  and  y  is  the  1 ine-of-sight  path 
of  a  source  ray  through  the  medium  (see  Appendix  I). 

For  a  highly  ionized  plasma  in  which  the  wave 
frequency  is  much  larger  than  the  electron  plasma  frequency, 
which  in  turn  is  much  larger  than  electron  cyclotron 
frequency  and  electron-ion  collision  frequency,  the  term  in 
the  integral  of  Equation  (1)  can  be  approximated  by 


M(y)-1 


-yy>2 
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(see  Appendix  II)  where  the  plasma  frequency  is  given  by 


The  fringe  shift  observed  on  the  screen  is  therefore  a  direct 
measure  of  plasma  density. 

In  this  research  the  critical  density  for  which  the 
plasma  frequency  is  equal  to  the  source  frequency  CUV  light, 
X=337. lnm)  is  9. 8‘ 102^cm-2,  which  represents  the  theoretical 
upper  limit  in  plasma  density  at  which  the  laser  light  is 
reflected.  In  practice,  the  upper  limit  is  about  one-third 
of  this  valued  From  Equations  (1),  (2)  and  (3),  the  minimum 
density  that  can  be  measured  in  this  experiment  is 
approximately  1017cm'2,  which  corresponds  to  a  few  tenths  of 


II.  EQUIPMENT 


The  interferometer  system  (Figure  3)  can  be 
separated  into  three  major  components: 

1.  a  radiating  source, 

2.  interferometric  optics,  and 

3.  a  detection  stage- 

A  laser  is  the  obvious  choice  for  the  radiating  source  since 
it  is  inherently  coherent  and  its  frequency  can  be  much 
greater  than  the  plasma  frequency  of  a  dense  plasma.  The 
interferometric  optics  consist  of  two  beamsplitters  (one  to 
split  the  light  beam  and  one  to  recombine  it)  and  two  beam 
guidance  mirrors.  A  special  camera  is  used  in  the  detection 
stage  to  record  the  fringes  on  film,  and  a  fast  photodiode  is 
used  to  determine  the  timing  of  the  laser  output  relative  to 
the  plasma  event. 


A.  The  Nitrogen  Laser 

A  molecular  nitrogen  laser  has  been  selected  as  the 
source.  The  wavelength  of  the  most  powerful  N9  line®  is 
X=337.  1  nm.  which  corresponds  to  the  C®II,,  -  BJIIC>  band. 

The  laser  was  built  at  the  University  of  Illinois 
and  was  originally  designed  to  be  used  in  an  undergraduate 
laboratory  as  a  learning  aid.  Several  modifications  were 
made  to  adapt  the  laser  to  the  interferometer,  such  as  a 
larger  electrode  cross  section  to  produce  a  more  uniform  beam 
and  a  low  Jitter  triggering  system. 


1.  Source  Requirements 

The  source  requirements  are: 

1)  the  frequency  must  be  high  enough  that  the 
plasma  acts  as  a  lossless  dielectric, 

2)  the  frequency  must  be  low  enough  to  retain 
sensitivity  in  the  interferometer  system, 

3)  the  output  must  have  a  fair  degree  of  coherence 
in  order  to  obtain  fringes, 

4)  the  source  should  emit  paraxial  rays, 

5)  the  source  must  have  a  high  intensity  output, 
and 

6)  the  source  should  be  pulsed  in  the  case  of  a 
dynamic  plasma. 

The  nitrogen  laser  fulfills  all  the  requirements  needed  for 
its  use  in  this  experiment.  The  first  requirement  is  dealt 
with  in  Appendix  II. 

The  second  requirement  arises  from  the  fact  that 
the  fringe  shift  is  inversely  proportional  to  the  source 
frequency  (Equations  (1)  and  (2)).  As  the  frequency 
increases,  the  fringe  shift  (for  a  given  plasma  density)  will 
decrease  and  the  sensitivity  of  the  system  will  be  reduced. 
The  useful  density  range  (see  Introduction)  is 

lO^crn'^  <  np  <  lO^cm-^  . 

The  maximum  density  for  the  University  of  Illinois  DPF  has 
been  estimated-*  to  be  ne-“2'  lO^crrf -b  Thus,  the  nitrogen 
laser  will  resolve  the  plasma  density  over  two  orders  of 


magnitude. 


The  third  requirement,  coherence,  is  needed  to 
obtain  fringes.  In  general,  the  optical  path  difference 
between  the  two  arms  of  the  interferometer  (including  the 
phase  shift  from  the  plasma  in  the  "probe"  path)  must  be  less 
than  the  coherence  length  of  the  source.  The  coherence 
length  of  a  nitrogen  laser  is  only  a  few  millimeters,  since 
the  laser  operates  without  a  cavity  (due  to  the  high  gain  of 
N9).  The  short  coherence  length  adds  inconvenience  to  the 
adjustment  of  the  path  length  difference  (for  further 
details,  see  Section  B:  The  Interferometer). 

The  fourth  requirement,  paraxial  rays,  is  needed  so 
the  phase  shift  is  not  lost  due  to  finite  detector  aperture. 
The  N2  laser  operates  without  a  cavity  and  the  beam 
divergence  is  relatively  large-  The  divergence  can  be 
approximated  by 

6/2  *  w/d  (4) 

where  w  is  the  cavity  width  and  d  is  the  electrode  length. 
For  the  N2  laser  used  in  this  research,  w=l.2  cm,  d=91  cm  and 
9~26  mrad-  A  pinhole  is  used  to  remove  most  of  the 
nonparaxial  rays  (see  Section  C:  The  Detection  Stage). 

The  intensity  of  the  source  (fifth  requirement) 
must  be  much  higher  than  the  intensity  irradiated  by  the 
plasma  at  that  wavelength-  The  light  from  the  plasma  adds  a 
back.grcund  component  to  the  fringe  visibility,  which  tends  to 


swamp  the  fringe  pattern  with  light.  The  output  power  of  the 
laser  is  ~500  kW  (Figure  4)  which  is  much  greater  than  the 
expected  blackbody  radiation  of  the  pinch  C~  watts). 


Finally,  the  laser  must  be  pulsed  since  the  plasma 
produced  in  a  DPF  is  rather  short  lived.  It  is  known13  from 
streak  pictures  taken  of  the  University  of  Illinois  DPF  that 
the  radial  implosion  velocity  of  the  pinch  at  a  fill  pressure 
of  3  torr  is  ~11  cm/ps.  Since  the  laser  pulse  is  6  ns  at 
FWHM  (Figure  4),  the  horizontal  resolution  of  the  fringe 
profile  is  ~0.7  mm.  The  resolution  is  actually  much  better 
because  a)  the  radial  velocity  represents  a  maximum3  and  b) 
the  fringe  smear  caused  by  the  plasma  movement  is  averaged 
out  when  the  data  are  read. 


Figure  4.  The  N?  laser  output  as  a  function  of  time 
(detected  using  a  fast  silicon  photodiode). 
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Laser  Operation 


Table 

1  lists  the 

specif icat ions 

of 

the  nitrogen 

laser. 

F igure 

5  shows  a 

block,  diagram 

of 

the  laser's 

electrical  system  and  Figure  6  is  of  the  gas  flow  system. 

The  energy  storage  system  for  a  nitrogen  laser  must 
be  able  to  excite  the  gas  to  the  C^IIU  state  in  a  time  much 
less  than  40  ns,  which  is  the  spontaneous  lifetime  of  the 
state.  To  achieve  this  it  is  common  to  use  a  Elumlein 
structure,  which  is  a  parallel-plate  transmission  line  with 
low  inductance.  As  shown  in  Figure  7a,  the  total  inductance 
and  capacitance  of  the  discharge  system  can  be  modeled  in  an 
LC  circuit  with  a  period  of 

T  =  2k(LC)^  (5) 

The  highest  potential  across  the  cavity  electrodes  is  reached 
at  1/4  of  this  period  (Figure  7b),  at  which  time  the  N9 


discharge 

takes  place. 

The  total  capacitance  on 

the 

osci 1 latory 

side  is  ~10  nF 

A  discharge  time 

of  much 

less 

than  40  ns  requires  a 

total  inductance 

of  L«  16 

nH. 

Therefore, 

the  design  of 

a  low  inductance 

spark  gap 

is 

crucial  since  this  is  probably  the  largest  source  of 
inductance  in  the  discharge  system.  A  field  distortion  gap 
with  a  "ring"  trigger  was  designed  to  fire  the  laser.  The 
size  of  the  gap  as  well  as  the  electrode  spacing  were  kept  to 
a  minimum.  The  gap  is  typically  pressurised  with  air  to 
~30  lb. 


Energy/Pulse : 
FWHM: 


Table  1.  Nitrogen  Laser  Specifications 


Population  Inversion  Method 


E/P: 


Electron  collisions  by 
transverse  electric 
amplification  CTEA) 

210  V/cm-torr 

~  3.0  mJ 

-  6  nS 


Energy  Storage: 


Ef  f iciency : 

Gas  Fill  System: 
Trigger  Type: 


Cavity  Type: 

Electrode  Separation: 


Electrode  Length: 


Cavity  Length; 


Double  Blumlein 
O  4. 1J  @  Vo=20  RV) 

0.  07% 


Static  fill 

Field  distortion  spark  gap 
with  ring  trigger 

Superradiant 

1.  2  cm 

91  cm 

100  cm 
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SOLENOID 
3-WAY  VALVE 


STATIC  FILL  LINE 


Figure  0.  End  view  of  No  laser  and  the 
gas  flow  arrangement. 


It  is  also  critical  that  the  jitter  of  the  laser 
pulse  with  respect  to  the  pinch  time  be  kept  minimal.  This 
was  accomplished  with  the  aid  of  a  thyratron  pulse  unit  used 
to  trigger  the  spark  gap.  As  seen  in  Figure  5,  the  thyratron 
and  the  Blumlein  are  charged  with  the  same  power  supply. 
However,  the  thyratron  releases  a  -20  kV  pulse  to  the  +20  kV 
charged  spark  gap  helping  to  reduce  Jitter  since  the 
breakdown  occurs  from  the  trigger  ring  to  both  main 
electrodes  simultaneously.  It  was  experimentally  determined 
that  a  small  inductance  of  %13  pH  is  needed  between  the 
thyratron  pulser  and  the  gap's  trigger  ring  to  reduce  the 
amount  of  energy  flow  to  the  thyratron  when  triggered. 

B.  The  Mach-Zehnder  Interferometer 

The  Mach-Zehnder  interferometer  consists  of  a 
source  beamsplitter  (BS^),  two  beam  guidance  mirrors  (M^  and 
M2)  and  a  beam  recombination  beamsplitter  (BS2)  shown  in 
detail  in  Figure  8  (setup  for  DPF  density  measurements).  All 
these  components  are  oriented  at  45°  angles  to  the  initial 
source  beam  with  adjustment  provisions  made  through  gimbal 
mounts  (CM).  Adjustment  for  the  equality  of  the  two  beam 
paths  is  provided  by  a  translational  stage  attached  to  mirror 
Mo-  This  also  eliminates  the  need  for  a  phase  shift 
compensator  in  the  reference  path  (to  account  for  the  phase 
change  due  to  the  windows  in  the  probe  path).  All  gimbal 
mounts  and  their  respective  optics  are  attached  to  the  DPF 


vacuum  chamber  with  flanges  designed  to  reduce  room  vibra¬ 
tions  that  may  be  present. 

Each  optical  component  that  the  beam  traverses  must 
be  flat  to  much  less  than  a  wavelength  of  the  source  in  order 
to  exhibit  a  straight,  flat  vacuum  fringe  pattern.  Since 
each  wave  traverses  different  optical  components,  the 
flatness  prevents  different  phase  distortions  in  each  path. 

The  substrate  for  all  optics  is  synthetic  fused 
silica.  This  material  has  exceptional  UV  transmission 
properties  (~90%  at  \=337  ran)  and  is  very  resistant  to 
radiation  darkening  from  UV,  x-rays,  and  neutrons-*-^,  all  of 
which  are  produced  in  a  dense  plasma  focus.  Three  tier 
baffles  have  been  placed  in  the  DPF  adjacent  the  fill  chamber 
windows  (Figure  8)  to  further  reduce  darkening  of  the  windows 
from  electrode  metal  vapor  deposits.  The  deposits  could 
destroy  the  flatness  of  the  windows.  Each  baffle,  as  well  as 
most  other  nonoptical  components  to  the  interferometer 
system,  have  been  blackened  to  reduce  reflections  for  safety 
purposes  and  to  reduce  other  possible  distortions  in  the 
fringe  pattern  induced  by  reflections. 

It  has  been  mentioned  that  the  difference  in  path 
lengths  between  the  reference  wave  and  the  the  probe  wave 
should  be  minimised  (to  much  less  than  the  coherence  length 
of  the  source!  in  order  to  obtain  highly  visible  fringes. 
This  is  especially  true  since  the  nitrogen  laser  does  not  use 
a  resonating  cavity.  In  such  a  case,  the  bandwidth  of  the 


.V. 


wave  imposes  a  limit  on  how  rapidly  the  phase  can  vary,  and 


this  can  be  related  to  the  coherence  length  by^ 


lc  =  C/Ai’ 


where  AV  is  the  bandwidth  and  c  is  the  speed  of  light.  For 


the  nitrogen  laser,  the  coherence  length  is  on  the  order  of  a 


few  millimeters.  However,  the  contrast  of  the  fringes  varies 


strongly  over  this  length.  The  variations  can  best  be 


described  in  terms  of  the  visibility  function^,  V,  defined  as 


Ima:<  ^min 


Imax  +  ^min 


where  Imax  and  Imin  are  the  constructive  and  destructive  wave 


intensities,  respectively.  It  is  a  function  of  the  optical 


path  difference,  t,  bandwidth  of  the  source,  AV,  and 


intensities  of  the  two  waves. 


The  visibility  function  for  the  system  used  in  this 


thesis  peaks  at  t=0  and  rapidly  approaches  zero.  The  loss  of 


visibility  in  the  fringe  pattern  as  t  increases  or  decreases 


is  caused  by  one  wave  becoming  more  random  in  phase  with 


respect  to  the  other  because  of  line  broadening  of  tht 


'hrcugh  experimentation,  the  path  difference  must  be 


less  than  r=20  pm  in  order  to  obtain  high  contrast  fringe 


(Figure  9).  The  procedure  used  to  align  the  interferometer 


m 
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within  this  limit  is  discussed  in  Appendix  III. 

The  visibility  function  is  greatest  when  the  two 
wave  intensities  are  equal.  Since  the  beamsplitters  used  for 
this  experiment  are  very  nearly  50/50.  the  intensities  are 
equal  up  to  the  point  of  recombinat ion.  However,  due  to  the 
antireflection  coating  on  the  back  of  BS2.  more  than  60%  of 
the  reference  beam  gets  transmitted  to  the  detector  whereas 
only  50%  of  the  probe  beam  reaches  it.  The  problem  is  easily 
solved  by  placing  an  80%  transmission  neutral  density  filter 
in  the  reference  path.  The  filter,  also  manufactured  from 
synthetic  fused  silica,  has  an  index  of  refraction  different 
from  the  surrounding  air,  which  affects  the  path  length 
difference  and  must  be  taken  into  account. 


'igure  0-  Vacuum  interf erogram  showing  highly  visible 
fringes  from  a  well  adjusted  interferometer. 
The  nonparallelism  between  the  top  and 
bottom  fringes  is  caused  by  skew  paths- 
The  problem  adds  only  minimal  distortions 
to  horizontal  measurements. 
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C.  The  Detection  Stage 

Two  detectors  are  needed  for  this  interferometer 
system:  one  for  timing  purposes  and  the  other  to  record  the 
fringe  pattern.  The  timing  detector  was  built  around  a  fast 
silicon  photodiode  (used  to  obtain  the  trace  in  Figure  4). 
The  fringe  recorder  is  Polaroid  film  which  the  light  reaches 
after  passing  through  a  series  of  optics  needed  to  correct 
distortions  and  add  magnif icat ion  to  the  fringe  profile. 

1.  The  Timing  Detector 

The  timing  detector  consists  of  a  fast  silicon 
photodiode  C 1  ns  risetime)  in  the  circuit  of  Figure  10.  The 
peak  of  the  laser  pulse  from  the  diode  indicates  the  time  at 
which  the  fringe  photograph  is  taken  relative  to  the  plasma 
pinch  event.  Because  of  the  high  sensitivity  of  the  diode 
Cl  A/W  at  X=337  nm) ,  several  neutral  density  filters  were 
needed  to  attenuate  the  laser  light.  A  bandpass  filter 
CX=337±6  nm)  at  the  diode  input  isolates  the  laser  light  from 
noise  produced  by  both  room  and  plasma  light. 

2.  The  Fringe  Detector 

The  fringe  detector  (Figure  11)  consists  of  three 

parts : 

a)  a  multipurpose  lens  and  spatial  low-pass 
filter  C  p i nno 1 e ) . 

b)  neutral  density  and  bandpass  filters,  and 

c)  photographic  film. 


BATTERY 
TEST  POINT 


0.  1  hF 


a. 


The  Lens  and  Pinhole 


The  most  important  purpose  of  the  lens  is  to 
prevent  distortion  of  the  fringe  image  from  laser  beam 
refraction  created  through  density  variations  of  the  plasma 
pinch.  Using  the  ray  matrix  method-1-1,  the  input  and  output 
parameters  of  a  general  optical  system  are  related  by  the 
ABCD  matrix, 


Rout 

A 

B 

Rin 

R  'out 

C 

D 

R  '  in 

where  R  is  the  radial  distance  from  the  optical  axis  and  R' 
is  the  slope  of  that  ray.  For  an  optical  system  consisting 
of  a  length  of  free  space  Cd^),  a  lens  of  focal  length  f 
fusing  thin  lens  approximations)  and  another  segment  of  free 
space  (d2),  the  matrix  becomes 


R2~ 

l-d2/f 

di~dj_d2/f  +d2 

1 

fH 

L* 

R<2 

-l/f 

l-dx/f 

LR'lJ 

By  letting 


C  9) 


B  =  d1-d1d2/f+d2  =0  C 10) 

where  is  the  distance  from  the  center  of  the  plasma  pinch 
to  the  lens  and  do  is  the  distance  from  the  lens  to  the  film 
plane,  the  radial  function  of  the  image  on  the  film  does  not 
depend  on  the  slopes  of  the  rays  as  they  leave  the  piasma. 
In  effect,  the  plane  of  the  film  is  moved  to  a  distance  B 
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from  the  plasma  pinch.  When  B=0,  the  system  is  in  focus. 

The  second  purpose  for  the  lens  is  to  magnify  the 
fringe  image  on  the  film  for  better  readability.  The 
magnification  factor  is  (for  B=0) 


!/  --  R2/^i  =  l-d2/f. 


C 11) 


By  choosing  a  particular  magnification  factor  and  focal 
length,  the  distances  and  d9  are  determined. 

The  third  purpose  of  the  lens  is  to  focus  the  near- 
paraxial  beam  through  a  sufficiently  small  pinhole.  The  term 
near-paraxial  is  used  in  reference  to  the  refraction  of  the 
laser  light  through  the  plasma  rather  than  from  the  natural 
beam  divergence  which  can  be  assumed  to  equal  zero  from 
Equation  CIO).  Since  Equation  CIO)  can  equal  zero  for  one 
image  plane  distance  only  Cd2=distance  from  the  lens  to  the 
film  plane),  a  particular  ray's  parameters  at  the  pinhole 
will  depend  on  its  angle  of  departure  from  the  plasma.  This 
applies  to  any  nonparaxial  ray  whether  it  be  from  the  plasma 
generated  light  or  the  laser  beam.  From  Equation  C9),  it  can 
be  seen  that  a  pinhole  distance  equal  to  the  lens  focal 
length  will  force  the  A  term  to  zero.  Thus,  the  nonparaxial 
rays  CR'^O),  most  of  which  are  from  the  plasma  light,  are 
blocKed  by  the  pinhole. 

In  order  to  predict  the  pinhole  size  needed  to  let 
the  refracted  laser  rays  through  but  still  block  most  of  the 
plasma  light,  do  in  Equation  C9)  is  defined  to  be  the 
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distance  from  the  lens  to  the  pinhole  and  the  maximum  R2  is 
calculated.  For  a  plasma  that  is  cy 1 indrical ly  symmetric  in 
its  density  profile,  and  R ' ^  can  be  related12  by 


0(b)  =  n-2  r  2[ b~2-r  2-b  2V(r)l~^  dr 


C  12) 


where 

R  =  tan  0(b) 

~  b  =  impact  parameter  of  the  given  ray 
0(b)  =  ray's  angle  of  departure  from  plasma 
V(r)  =  tfpCr)/fwave]2 

and  rQ  is  the  ray's  distance  of  closest  approach  to  the 
plasma  center  (calculated  by  setting  the  expression  in 
brackets  in  Equation  (12)  equal  to  zero).  Shmoys12  presents 
a  parabolic  distribution  for  V(r): 


V(r)  =  VQ [ l-( r/a)^l ,  r<a 
V(r)  =  0,  r>a 


The  resulting  analytical  solution  for  0(b)  is 


(13) 


VQ  [  (b/a)2  -  (b/a)4] 

0(b)  =  sin-1 - — — 

{V0(b/a)2+(( l-V0)/2] 2V* 


(  14) 
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For  a  peak  density  of  ne  =  2-  10igcm-2,  VQ  ~  0.002.  Also, 
from  streak  photographs  of  a  typical  pinch  produced  by  the 
University  of  Illinois  DPF,  the  "a"  parameter  in  Equation 
(13)  is  estimated  to  be  approximately  1-5  mm.  Using  these 
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values  in  Equation  C14)  and  then  in  Equation  (9),  a  pinhole 
size  of  2  mm  diameter  is  estimated.  However,  the  2  mm 
pinhole  also  blocks  much  of  the  diverging  light  from  the 
laser  light.  This  in  turn  reduces  the  picture  area  of  the 
pinch  so  a  pinhole  of  3  mm  diameter  was  selected  to  increase 
the  picture  size  with  little  effect  on  fringe  visibility. 

Since  a  simple,  single  lens  system  is  used,  it  is 
important  to  reduce  as  many  lens  aberrations  as  possible  by 
lens  shape.  Due  to  the  nature  of  the  laser  light,  only 
monochromatic  aberrations  are  important.  Those  most  affected 
by  lens  shape  are  spherical  aberration  and  coma.  The  former 
type  is  the  failure  of  rays  from  a  point  object  on  the 
optical  axis  to  converge  to  a  point  image.  They  actually 
converge  to  a  small  circle  at  the  image  plane  which  indicates 
the  finest  detail  that  the  lens  can  resolve.  Coma  is  a 
similar  effect  but  for  off-axis  rays  which  converge  to  a 
comet-shaped  configuration.  For  the  case  of  paraxial  rays 
merging  to  the  focal  point  of  the  lens,  these  aberrations  can 
be  minimized  by  choosing  a  plano-convex  lens  with  the  convex 
side  toward  the  paraxial  light  ray  source^. 

The  shape  and  index  of  refraction  of  a  lens  affect 
the  focal  spot  size  through  a  function  ECp)  where  M  is  the 
index  of  refraction  for  a  given  wavelength  of  light.  From 
third-order,  thin- lens  aberration  theory  and  from  the 
Fraunhofer  diffraction  limit  for  a  uniformly  illuminated 
square  aperture  of  side  dimension  a,  the  spot  size  diameter 


can  be  approximated  by 


s  *  2fCK(p)/2  (D/f)J  +  X/al  (15) 

where  D  is  the  larger  dimension  and  a  is  the  smaller 
dimension  of  the  beam  at  the  lens.  For  a  plano-convex  lens*® 


K.(M)  = 


C32CM- l)2  1  • 


From  Equation  (15)  the  optimum  lens  f-number  (f/D)  can  be 
found  for  which  the  spot  size  is  minimized: 


(  15) 
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Cf/D)opt  =  t3fKCfO/2X]*. 


C  17) 


For  D  =  30  mm,  X  =  337. 1  nm  and  M  =  1.47977  (optical  quality 
synthetic  fused  silica),  we  have  a  focal  length  f=616  mm  and 
a  spot  size  s=22  Pm.  Due  to  the  smaller  cost  of  standard 
items,  a  lens  of  f  =  500  mm  and  diameter  of  50  mm  was  chosen. 
As  with  the  other  optics  in  the  system,  high  quality  in 
manufacturing  is  needed  to  reduce  irregular  aberrations  such 
as  scratches,  digs,  surface  contamination  and  inhomoge¬ 
neities. 

For  this  particular  interferometer,  a  magnif icat ion 
factor  of  M=2  was  chosen  due  to  physical  constraints. 
Therefore,  from  Equations  (10)  and  (11),  the  distance  from 
the  plasma  pinch  to  the  lens  must  be  750  mm,  the  distance 
from  the  lens  to  the  oinhole  must  be  500  mm,  and  the  distance 
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magnification  factor  is  M=2.  2  both  horizontally  and 
vertically  (Figure  12).  The  extra  0-2  in  M  is  within  the  lens 
focal  length  tolerance  (±5%). 


Figure  12- 


Vacuum  interf erogram  with  the  image  of  a 
washer  placed  at  the  center  electrode  used 
to  check  focus  and  hor./vert.  magnification 
(Mh>  v  =  2.  2) . 


b-  Filters 


The  second  component  in  the  fringe  detection  stage 
is  the  optical  filter  system.  A  combination  of  neutral 
density  filters,  each  of  a  different  transmittance,  was  used 
to  experimentally  match  the  laser  light  intensity  to  the  film 
sensitivity  in  order  to  obtain  the  highest  visibility  of  the 
fringe  pattern.  Also  a  337.1  nm  bandpass  filter  of  bandwidth 
-6  nm  is  placed  before  the  film  to  reduce  plasma  and  room 
light.  This  filter  should  be  placed  as  close  to  the  film  as 
possible  in  order  to  preserve  the  phase  information  within 
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the  beam.  No  such  problems  occurred  with  the  particular 
neutral  density  filters  used.  These  filters  were  placed  near 
the  pinhole.  A  dark  canopy  was  placed  over  the  detector  from 
the  pinhole  to  the  camera  back  covering  the  intermediate 
filter  system  to  prevent  room  light  exposure. 

c.  Film 

Originally,  the  film  chosen  to  record  the  fringe 
pattern  was  Polaroid  type  665-negative  pack  film.  Pack  film 
was  selected  to  assist  in  faster  DPF  data  taking  and  for 
general  convenience.  Type  665  was  selected  because  of  the 
slow  speed  and  therefore  smaller  sensitivity  to  background 
room  light.  Also,  the  photograph  negatives  have  much  higher 
resolution  than  their  associated  positives  C 170  line  pairs/mm 
as  compared  to  17  line  pairs/rnm  for  its  positive 
counterpart).  A  problem  with  the  negatives  is  that  they 
darken  even  in  dim  room  lighting,  which  hampers  the  fringe 
readability.  Also,  the  film  is  relatively  expensive  and  hard 
to  obtain.  Because  of  these  reasons,  Polaroid  type 
667-positive  pack  film  was  used.  This  film  is  comparable  to 
the  665-positive  in  resolution  but  is  *40  times  faster. 
Neutral  density  filters  with  lower  transmittance  were  used  to 


attenuate  the  beam. 

D.  The  Complete  Interferometer  System 

The  complete  interferometer  system  is  shown  in 
Figure  13.  For  ease  in  construction,  the  laser, 
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interferometer  optics  and  detectors  are  not  physically 
connected.  The  laser  and  fringe  detector  rests  on  a  brick 
wall  with  a  1/4  inch  rubber  mat  insulating  the  components 
from  possible  floor  vibrations  Csuch  as  from  mechanical 
pumps)  preventing  the  need  for  an  expensive  optics  table. 
The  large  mass  of  the  wall  aiso  protects  the  laser  and  fringe 
detector  from  most  air  vibrations  that  could  disrupt  critical 
alignment  of  the  system. 
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III.  DATA  ANALYSIS 


Once  the  fringe  profile  of  the  plasma  is  obtained, 
the  density  can  be  found  through  Equations  Cl).  (2)  and  (3). 
Symmetry  conditions  are  assumed  to  simplify  the  analysis. 
The  raw  data  are  numerically  transformed  into  a  density 
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Analytical  examples  are  given  as  a  test  of  the 


software  and  for  numerical  error  analysis. 


A-  Line  Density 

Line  density  is  the  number  of  particles  per  unit 
length.  This  parameter  is  useful  for  the  study  of  particle 
conservation  anc,  axial  movement  in  a  dynamic  plasma. 

The  optical  phase  shift  of  a  ray  traversing  a  chord 
through  a  plasma  (with  respect  to  vacuum)  is  given  by 


CIS) 


S (x)dx  =  —  [p(x,y)-l]  dy 


where  £(x)  is  in  wavelength  units.  This  is  the  same  as 
Equation  Cl)  (the  x  dependence  is  included  here).  Through 


the  us1; 


of  Eouations  (2)  and  (3),  Equation  (18)  can  be 


reorganised  to  yield 


neCx,  y )dy  =  -DXrySCx) 


(19) 


where  nc  is  the  critical  density  where  the  plasma  frequency 
is  equal  to  the  wave  frequency.  The  integration  of  both 
sides  over  the  x-axis  will  directly  give  the  number  of 
particles  per  unit  length  along  the  z-axis  of  the  plasma. 
Therefore,  the  line  density  (n^)  can  easily  be  calculated 
from  the  fringe  shift  data  by  using 


n^  =  2Xnc 


S (x)  dx 


B.  The  Plasma  Density  (Assuming  Axial  Symmetry) 

Given  the  geometry  of  the  DPF,  a  cy  1  indrical  ly 
symmetric  plasma  is  assumed  (r2=x2+y2)-  For  a  ray  of  width 
dx  traversing  a  plasma  slice  of  height  dz  with  an  impact 
parameter,  x=b  (Figure  14),  Equation  (19)  becomes 


£(x)  - 


n(  r ) 

■J  T2 


-r  dr 


where  R  is  the  outer  boundary  of  the  plasma  (density  is 
negligible  for  all  r5:R).  Equation  (21)  is  a  form  of  the  Abel 
integral  and  can  be  inverted  to  yield 


ne(r)  - 


d£(x)  dx 
dx  -iyj  -v‘ 
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Figure  14. 


Relationship  between  fringe  pattern 
and  top  profile  of  cy 1 indrical ly 
symmetric  plasma  used  for  numerical 
analysis.  R=aN;  n=2q. 
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This  is  the  form  which  is  needed  to  obtain  the  radial  density 
profile  of  the  plasma  from  the  fringe  shift  data  on  the  film. 


C.  Numerical  Analysis 

In  general,  numerical  differentiation  is  much  less 
accurate  than  numerical  integration14.  The  data  from  the 
fringe  pattern  are  read  in  tabular  form.  As  seen  from 
Equation  (22),  the  data  must  be  differentiated  and  any  noise 
would  be  greatly  amplified,  especially  at  large  gradients  in 
the  phase  shift  curve. 

William  L.  Barr15  introduced  a  method  to  numeri¬ 
cally  solve  equations  similar  to  Equation  (22).  In  his 
method.  the  order  of  differentiation  and  integration  is 
reversed,  and  a  least-squares  smoothing  technique  is  applied 
before  the  final  differentiation.  The  only  assumptions  used 
in  his  method  are  that  the  curve  to  be  inverted  is 
axisymmetric  and  that  its  slope  on  axis  is  zero. 

The  accuracy  of  Barr 's  method  depends  on  the  number 
of  data  points  used.  The  use  of  twenty  data  points  results 
in  good  accuracy  for  Gaussian  and  polynomial  distributions1^ 
Such  profiles  have  been  observed  in  other  plasma  pinch 


experiment 


.17,  18.  19 


Barr  divides  the  x-axis  into  equal  increments  a, 


such  that 


xn  =  na,  rk  =  ka.  R  =  Na 


(23) 
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where  n  and  k  are  integers  between  0  and  N,  inclusively.  The 
data,  however,  are  most  easily  read  in  equal  increments  of 


S(x).  In  order  to  obtain  the  required  SCx^,  a  third-degree 
Lagrangian  interpolating  polynomial14  is  used  where 


cf(xn)  =  *C:<o> 


+  6  a x1) 


Cxn-x1)(xn-x2)Cxn-x3) 
(Xq-X1) (Xq-X?) (Xq-Xj) 

(xn‘x0)(xn'x2)(xn-x3) 
(X-^-Xq)  (x1-x2)  (X1-X3  ) 


Cxn-x0)Cxn-Xi)Cxn-x3) 

+  6 ( Xo ) - 

^  CX2-Xq) Cx2-x1) (X2-X3) 


C  24) 


+  <SCx3) 


Cxn-Xo)CXn-x1)Cxn-x2) 

(X3-Xo)Cx3-x1)(x3-x2) 


Figure  15  demonstrates  the  method  with  a  Gaussian  distri¬ 
bution.  A  Gaussian  was  selected  since  a  single  polynomial 
cannot  exactly  fit  it.  The  interpolation  error  (excluding 
roundoff  error)  is  within  0.2%  of  the  peak  in  the  curve. 

For  a  near  symmetric  £(xn)  distribution,  the  points 
on  either  side  of  the  pinch  are  averaged  in  order  to  obtain 
the  required  zero  slope  at  xn=0.  If  the  data  are  not  sym¬ 
metrical,  the  density  on  either  side  of  the  pinch  axis  can  be 
separately  predicted  by  mirror  imaging  the  data  points  of  one 
side  over  to  the  other  side  of  xn=0.  If  the  data  are  sym¬ 
metric  but  off  the  chosen  pinch  axis  (such  is  the  case  for  a 
m=l  kink  instability),  the  axis  of  symmetry  is  taken  as  x=0 
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Figure  15. 


Ca)  Diagram  showing  the  transformation 
of  the  fringe  shift  data  to  equal  increments 
in  (b)  Input  data  to  the  Lagrangian 
interpolating  polynomial  of  equal  increments 
in  S ( x)  • 


mmam mmmsm 


40 


for  the  analysis. 


Once  the  £(xn)  values  are  determined,  Barr's  method 


of  inverting  the  data  will  smooth  any  round-off  errors.  The 
final  form  of  Equation  (22)  then  becomes 


neCxk)  = 


2 Xr^  N 

1  (W^n^ 


n=k-2 


k  >  2 


(25) 


2  An,-.  N 

ne(xk)  =  -  2  8i/n'  £ (xn) , 

e  K  n=0  -n  n 


k  <  2 


where  0kn  are  constants  that  include  the  entire  integration, 
least-squares  fitting,  and  the  final  differentiation  process 


of  the  <S(xn)  curve.  The  values  are  listed  in  Table  2. 


D.  Analytical  Examples  and  Numerical  Error  Analysis 

The  equations  from  the  last  section  are  included  in 
a  computer  program.  The  software  was  written  to  accept  the 
raw  data  read  from  the  photographs  as  input.  The  phase  arid 
density  profiles  are  calculated  and  graphically  displayed. 

The  input  data  are  expected  to  be  between  <S(-x)=0 
and  <^(x)=0  with  the  axis  of  the  pinch  at  x=0.  Although  these 
endpoints  are  not  used  in  Barr's  method  of  inversion,  they 
are  included  in  the  Lagrangian  interpolating  polynomials  for 
better  accuracy  of  <S(x)  near  the  endpoints.  The  data  are 
interpolated  to  obtain  20  equidistant  points  on  each  side  of 
the  axis  where  the  separation  is  given  by  A=xmax/21  where 


xmax  is  the  larger  absolute  x  value  from  £(-x)  and  £(x). 


This  is  a  prerequisite  of  using  Barr's  method  of  inversion 
(here,  R=20A).  The  interpolated  data  on  either  side  of  the 
phase  profile  axis  (x=0)  are  then  averaged  in  order  to  obtain 
the  cy 1 indrical ly  symmetric  profile  in  r  which  is  required 
for  Abel  inversion.  After  the  inversion  is  performed,  the 
density  magnitude  at  each  point  r^=kA  (0<k<20)  is  acquired. 

To  obtain  a  more  uniform  display  of  the  density 
graphics,  the  inverted  data  are  again  interpolated,  using 
third-degree  Lagrangian  polynomials  to  fill  in  the  spaces 
between  r^.  The  program  is  straightforward  and  is  not  listed 
here. 

To  test  the  above  method,  two  representative  func¬ 
tions  are  assumed,  and  their  analytical  solutions  are 
compared  to  the  program  results.  The  first  function  is  a 
Gaussian  phase  profile-- 


ffCx)  =  e"x2/a2  .  C  26) 

A  Gaussian  profile  is  a  good  test  since  a  low-order  poly¬ 
nomial  cannot  exactly  fit  it  and  because  it  is  sharply 
peaked.  The  analytical  solution  for  the  Abel  inversion  of 
Equation  C26)  Cas  R  approaches  <*>)  is 
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Although  the  numerical  method  can  not  approach  this  limit 
(£(r)=0  for  r>R),  the  error  is  negligible  for  large  enough 
values  of  R.  The  numerical  and  analytical  inversions  of 
Equation  C26)  are  compared  in  Figure  16  for  a-V2.  The  maximum 
error  is  within  2%  of  the  peak  value  of  which  1%  is  from  the 
round-off  to  2  decimal  places. 

The  second  test  function  is  a  third-order  poly¬ 
nomial  in  y?- .  This  curve  represents  a  hollow  density  pro¬ 
file-  The  function  can  be  written,  in  general,  as 


(28) 


SCx)  =  a+bx^-cx^+dx^ 


The  solution  for  the  analytical  inversion  is 


A(r)  =  -  -vR^-r^  [  C  24d)  r4  +  (  12dR2-20c )  r2  +  C  9dR4-  10cR2  +  l5b)  1  . 

15tt 

(29) 


The  constants  in  Equation  (28)  are  chosen  such  that  the 
magnitude  and  slope  of  <S(x)  at  x=l  are  equal  to  zero.  This 
makes  the  obvious  choice  for  R=l.  For  a=l,  b=2,  c-7  and  d=4, 
Equation  (29)  becomes 


ACr)  =  -  41-r*  [ 24r4-23r“- 1 1 

ltm 


(30) 


In  Figure  17,  the  numerical  method  is  compared  to  Eq.  (30). 
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(a)  Comparison  between  the  numerical  and 
the  analytical  Abel  inversions  for  a 
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Figure  17. 
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(a)  comparison  between  the  numerical  and 
the  analytical  Abel  inversions  for  a 
third  order  polynomial  in  x;  (max.  error 
=  4%).  (b)  Program  input  data,  incremental 

in  6 (x) • 
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The  large  error  in  the  numerical  curve  near  r=l  is  caused  by 
the  sharp  curvature  and  low  values  in  £(x)  in  this  region. 
However,  the  curve  generated  numerically  is  meaningless 
beyond  r=0. 95  since  A(r)  does  not  go  to  zero  at  r=1.0.  The 
corresponding  points  are  not  used  in  the  experimental  data 
inversions.  The  maximum  error  generated  through  the 
numerical  inversion  is  therefore  within  4%  of  the  maximum 
value  in  the  region  of  interest. 
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IV.  DATA 

The  plasma  formed  by  the  dense  plasma  focus  is 
short  lived  (typically  <  100  ns).  The  jitter  of  the  laser 
trigger  must  be  minimized  in  order  to  accurately  determine 
the  timing  of  the  interferogram  exposure  with  respect  to  the 
dynamic  plasma.  The  last  section  in  this  chapter  is  on  error 
margins  for  the  measured  densities. 

A.  Timing 

The  jitter  time  of  the  pulse  used  to  trigger  the 
oscilloscope,  and  therefore  the  laser  (see  Figure  5  in 
Chapter  II),  is  critical  since  the  plasma  changes  rapidly 
with  time.  The  best  start  pulse  signal  seems  to  be  given  by 
the  voltage  induced  in  a  small  magnetic  probe  placed  between 
the  inner  and  outer  electrodes  (also  used  to  determine  the 
pinch  current).  The  trigger  system  then  bypasses  any  Jitter 
created  during  the  initial  breakdown  of  the  gas  along  the 
base  insulator  of  the  DPF  but  is  also  early  enough  to 
accommodate  for  electronic  equipment  delay. 
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recent  comparison  of  x-ray  peaks  with  laser  pulses  confirms 
the  timing  reported  in  the  data  section. 
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The  dynamic  behavior  of  the  plasma  in  the  DPF  can 
be  separated  into  three  stages  (.Figure  18): 


1.  first  compression. 

2.  radial  expansion  and 


3.  second  comoression  (olasma  breakuo). 


The  first  compression  lasts  for  approx  i  mat  iy  15  ns  in  the 
University  of  Illinois  dense  plasma  focus  at  the  end  of  which 
a  dense  plasma  column  is  formed.  The  column  then  bounces 
radially  outward  (radial  expansion  phase).  This  phase  can 
last  for  ^20  ns.  Finally,  a  second  compression  of  the  plasma 
occurs  in  which  the  pinch  column  breaks  up.  The  plasma 
breakup,  typically  caused  frc:..  .,.=0  (sausage)  instabilities, 
lasts  typically  less  than  25  ns. 

Included  in  this  section  are  representative 
interf erograms  and  their  respective  density  profiles  as  a 
function  of  r  and  a  (radial  and  axial  dimensions)  for  the 
three  pinch  phases  (although  from  different  shots).  The 
initial  fill  gas  for  all  data  is  deuterium  and  the  input 
energy  is  12.  5  kJ.  The  error  is  within  10%  of  the  maximum 
density  for  each  profile  along  the  s-axis  (see  Section  C.  of 
Chaoter  IV). 


The  First  Compression 

The  collapse  phase  is  recognised  by  the 
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RADIAL  EXPANSION 
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Cb) 


2nd  COMPRESSION 
CPINCH  BREAKUP ) 


Figure  13. 


stage 


(c) 


Phase  shift  data  of  (a)  the  first  compression 
of  the  pinch  cycle:  P=3.5  torr  (Ho).  V0  =  26  kV , 

(b)  the  radial  expansion  phase:  P=3. 2  torr  (Do). 

V.,  =  25  kV  and  (c)  the  second  compression  from  which 
the  plasma  is  breaking  up  due  to  m=0  instabilities: 
P=3. 5  torr  (Ho).  VQ=26  kV.  The  flat  at  the  bottom 
of  the  fringes  is  the  center  electrode.  It  is  noted 
here  that  sometimes  the  fringes  appear  inverted. 

This  is  due  to  the  sign  of  the  angle  between  the 
reference  and  probe  waves  upon  recombinat ion  and 
is  not  a  problem  for  analysis.  The  magnification 
for  all  interf erograms  is  2.2  unless  otherwise 
specif ied. 
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characteristic  shape  and  smooth  edges  of  the  fringe  profiles 
as  shown  in  Figures  19-21.  Although  these  edges  are  produced 
by  the  current  sheath  surrounding  the  plasma,  no  evidence  of 
dangerous  instabilities  is  present.  Figures  19  and  20  are  of 
the  early  collapse  phase.  Figure  19  is  earlier  than  Figure 
20  as  noted  from  the  larger  radius,  lower  density  and  hollow 
profile.  The  plasma  in  Figure  19  seems  to  be,  for  the  most 
part,  from  the  current  sheath,  although  some  new  plasma  seems 
to  be  forming  at  the  sheath  front. 

The  density  profile  of  Figure  20  suggests  that 
pinching  of  the  plasma  has  begun.  The  peak  density  is  more 
than  a  factor  of  three  greater  than  that  of  Figure  19.  Also, 
the  density  profile  is  no  longer  hollow.  The  flat  portion  of 
the  profiles  is  closed  to  that  of  the  uncompressed  plasma  of 
Figure  19  suggesting  this  portion  of  the  plasma  is  still 
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Both  figures  are  from  a  fill  pressure  of 


Figure  21  shows  the  end  of  the  first  compression 
when  the  oensity  is  maximum  and  the  radius  is  minimum.  The 
minimum  radius  occurs  between  2  mm  and  5  mm  from  the  hollow 
center  electrode  (3.5  mm  on  the  average;  and  is 
character ist ical ly  ^1  mm  ( FWHM )  as  previously  observed  with 
streak  photographs^-  of  the  visible  light.  The  shape  of  the 
profile  at  "  =  1.87  mm  in  Figure  21  and  the  high  density 
C'-IO1 -'em  suggest  that  the  olasma  is  at  maximum 


at  maximum 


compression. 


PLASMA  DENSITY 
C SCALE :  1 . 0= 1 .  6‘  1018CM-3 ) 
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(a )  Inter f erogram  and  cb)  electron  pi 
density  of  the  current  sheath  enteric 
first  compression  phase:  P=3. 2  torr  C 
V0=25  kV. 
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Figure  21. 

(a)  Interf erogram  and  (b)  electron  density 
profile  of  the  pinch  at  maximum  compression 
P=3. 9  torr  (Do),  Vn  =  25  KV. 
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Due  to  the  lack  of  x-ray  timing  for  the 
interf erogram  in  Figure  22,  and  since  this  plasma 
configuration  has  not  been  observed  since,  it  is  not  clear 
what  stage  the  pinch  was  in.  The  high  density  OlO^cm-3) 
and  small  minimum  radius  (<1  mm  at  2=2. 82  mm)  suggest  that 
the  pinch  is  near  maximum  compression. 

Disturbances,  similar  to  the  one  at  the  outer  left 
side  in  Figure  22,  typically  occur  in  pinches  before  maximum 
compression  (from  more  resent  x-ray  timed  interf erograms) . 
Disturbances  are  also  evident  at  the  right  side  of  the 
interferogram  of  Figure  21.  These  disturbances  do  not  occur 
at  the  minimum  radius  of  the  pinch  as  would  be  expected 
according  to  the  Kruskal-Schwarzschild  Theory23.  Such 
occurrences  have  been  reported  in  other  DPF/N2  interferometry 
experiments17,  18,  19  but  not  much  has  been  said  about  them 
because  they  disappear  during  the  plasma  expansion.  The 
disturbances  seem  to  be  weak  instabilities  of  long 
wavelength,  since  it  has  been  reported18  that  they  only 
appear  in  longer  plasma  columns.  From  analysis  of  the 
fringes  in  these  disturbances  Cnot  included  here),  the 
density  in  the  region  is  typically  of  low  magnitude  (^15%  of 
ne-max  in  F1§ure  22). 

The  average  line  density  over  several  first 
compression  data  is  ~7.  8'  1017/cm  at  the  tip  of  the  center 
electrode  and  ~1.4- 1018/cm  at  1  cm  above  the  center  electrode 
indicating  mass  flow  along  the  axis. 
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Figure  22.  (a)  Interf erogram  showing  a  disturbance  at 

the  outer  left  side  and  (b)  electron  density 
profile  of  the  plasma  thought  to  be  near 
maximum  compression:  P  =  3.  2  torr  (Do). 

V - 2 5  kV. 


SR 

i 


m 


r«w> 


i«*  i.«  i.<‘t.»  i.«'l 


2.  The  Radial  Expansion 

The  profile  of  Figure  23  at  2=1.41  mm  (minimum 
radius)  shows  a  wave  of  plasma  starting  its  radial  expansion 
phase  (the  "hump"  in  the  profile  resembles  the  profiles  in 
the  expansion  phase  of  Figure  24).  The  expansion  appears  to 
be  occurring  before  the  maximum  compression  at  larger  z 
(zippering).  Similar  pinch  profiles  have  been  used  in  the 
modeling  of  pinch  dynamics^. 

Figure  24  shows  the  plasma  at  maximum  expansion. 
The  radius  of  the  plasma  at  this  stage  is  typically  three 
times  the  minimum  radius  of  the  first  compression  and  the 
density  profile  is  hollow.  Also,  the  plasma  configuration  is 
generally  cylindrical  as  compared  to  the  curved  nature  of  the 
first  compression  stage.  The  edges  of  the  fringe  pattern 
tend  to  wiggle  rapidly  along  the  z-direction  indicating  short 
wavelength  instabilities  in  their  initial  stage1. 

Increasing  densities  along  the  pinch  axis  in  the 
profiles  of  Figure  24  also  suggest  axial  mass  flow  away  from 
the  center  electrode.  The  fill  pressures  for  the  expansion 
data  are  as  given  in  the  figure  captions. 


3.  The  Second  Compression  (Plasma  Breakup) 

The  fringe  data  for  this  stage  are  generally  not 
symmetr leal  and  cannot  be  inverted  to  get  the  density 
profiles.  Thus,  only  the  interf erograms  are  shown.  However, 
the  peak  densities  can  be  estimated  by  using  selected 
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The  instabilities  of  Figure  25  seem  to  have  grown 


to  the  dangerous  stage  just  before  pinch  breakup  occurs.  The 
wavelengths  of  the  instabilities  are  longer  than  the  rapid 
wiggles  of  Figure  24. 


Figure  25. 


Interf erogram  of  the  pinch  at  the  start 
of  the  second  compression.  The  m=0 
instabilities  responsible  for  the  plasma 
destruction  are  forming  at  the  edges  of 
the  fringes. 


Figure  26  consists  of  four  fringe  patterns  of  the  plasma  well 
into  the  breakup  phase.  These  photos  and  others  not 
reproduced  in  this  thesis  reveal  that  m=0  instabilities  are 
almost  always  responsible  for  the  plasma  destruction, 
although  m=l  (kink)  instabilities  have  been  observed.  From 
the  interf  erograrns,  the  wavelength  of  the  most  rapidly 
growing  instabilities  was  measured  and  varies  from  2  mm  to 
4  mm. 
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Figure  26. 


fa),  Cb).  (c)  and  (d)  are  four  inter- 
ferograms  showing  the  m=0  instabilities  that 
break  up  the  plasma  during  the  second 
compression.  P Q - J . 2  tori  and  V0-25  kV  for 
< a)  through  (d) . 


The  densities  in  the  constricted  regions  of  the 
instabilities  approach  2‘101^cm_^.  The  density  in  the 
"blobs"  is  around  4-  lO^cm--5,  and  the  line  density  is 
n]^10l7/cm,  which  is  an  order  of  magnitude  larger  than  in  the 
constricted  regions.  The  use  of  line  densities  also 
demonstrates  that  particle  conservation  is  not  preserved 
along  the  pinch  axis  between  stages.  It  is  believed  that  the 
particle  loss  is  mostly  out  the  ends  of  the  pinch  because  of 
the  z-component  of  the  kinetic  pressure  (away  from  rmin)  and 
lack  of  confinement  in  this  direction. 

C.  Error  Margins 

The  largest  error  incorporated  into  the  density 
profiles  comes  from  the  reading  of  the  fringe  patterns.  The 
data  were  read  in  1/2-fringe  increments  as  a  function  of  x 
using  a  vernier  caliper  of  0.1  mm  resolution.  By 
progressively  shifting  the  original  data  for  several  profiles 
from  x=0±0  mm  to  x=R±0.  1  mm  in  a  linear  fashion,  a  typical 
error  margin  is  calculated  to  be  less  than  ±3%  of  ne_max- 
The  method  was  chosen  since  a  particular  point  of  the  density 
profile  depends  only  on  the  larger  radius  points  in  the 
fringe  profile  when  the  data  are  Abel  inverted.  With  the 
addition  of  the  numerical  error  calculated  in  Chapter  III.  it 
can  be  argued  that  the  measured  densities  are  well  within 
±10%  of  n0_ma>..  From  Figjres  16  and  17  c;  Chapter  III,  it  is 
seen  that  the  measured  densities  underestimate  the  real 


values  at  the  peaks  and  overestimate  the  real  values  near  the 
edges.  The  error  in  the  line  density  calculations  can  be 
assumed  to  be  within  ±5%  of  the  real  values.  This  is  because 
the  Abel  inversion  was  not  needed  in  their  determination,  and 
the  same  polynomials  used  to  interpolate  the  data  were  also 
used  to  integrate  them. 
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V.  SUMMARY  AND  CONCLUSIONS 

The  purpose  of  this  work  was  to  measure  the 
electron  density  of  a  dense  plasma  focus  with  the  use  of 
laser  interferometry.  A  nitrogen  laser  was  used  as  the  probe 
beam,  since  it  fulfilled  all  of  the  requirements  imposed  by 
the  high  density  and  short  duration  of  the  plasma.  The  rest 
of  the  equipment  needed  for  this  research  has  been  discussed 
in  detail.  Also  discussed  were  the  numerical  methods  and 
assumptions  used  to  obtain  the  density  profiles  as  a  function 
of  the  radial  and  axial  dimensions. 

The  electron  density  in  the  University  of  Illinois 
plasma  focus  is  on  the  order  of  1019cm~3  at  maximum 
compression.  By  observing  the  different  phases  in  the  plasma 
event  leading  to  and  following  the  maximum  compression,  the 
dynamics  of  the  pinch  is  better  understood.  The  radius  of 
the  pinch  at  maximum  compression  is  ~1  mm  and  then  expands  to 
a  cylindrical  shape  with  radius  ~3  mm  before  plasma  breakup 
when  the  current  stops  flowing.  The  breakup  is  almost  always 
caused  by  m=0  instabilities  although  m=i  instabilities  have 
been  observed. 

The  Mach-Zehnder  interferometer  is  now  being  used 
to  understand  the  pinch  dynamics  when  impurities  are  added  to 
the  initial  gas  fill  (such  as  neon)  and/or  an  axial  magnetic 
field  is  added  at  the  pinch  axis-  Prel-.xnary  results  of  5% 
Me  added  to  ~3  torr  hydrogen  show  a  much  tighter  and 
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straighter  pinch  at  first  compression.  The  magnetic  field 
<~500  Gauss)  stabilizes  the  pinch  from  m=0  instabilities  but 
not  from  m=l  instabilities.  Perhaps  a  stronger  axial  field 
will  reduce  the  m=l  instabilities  and  prolong  the  confinement 
time  of  the  plasma. 
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APPENDIX  I.  FRINGE  SHIFT  OF  A  WAVE  IN  A  DIELECTRIC 

The  purpose  of  this  Appendix  is  to  demonstrate 
simply  the  relationship  between  fringe  shift,  £(y),  and  index 
of  refraction,  pCy),  of  a  lossless  dielectric  of  length  L  in 
one  arm  of  a  Mach-Zehnder  interferometer  (.'Figure  1.1).  The 
fringe  shift.  £(y),  is  found  for  one  ray  of  cross-sectional 
area  dx- ds  in  an  infinite  plane  wave. 

Written  in  phasor  notation,  a  wave  travelling  in 
the  y  direction  can  be  expressed  by 

E  =  Eq  e-J  C“t-0y)  (It 1) 

where  E0  is  the  constant  wave  amplitude,  w  is  the  radial  wave 
frequency  and  j3  is  the  phase  coefficient.  Equation  CI.l)  can 
similarly  be  written  as 

E  =  E0  e~J  Cvt-y>  Cl. 2) 

where  v  is  the  phase  velocity  of  the  wave.  In  fact, 

w  c 

v  =  -  =  - - -  C I.  3) 


medium.  The  index  of  refraction  is  related  to  the  relative 
permittivity  by 


Cl.  4) 


It  is  obvious  that  a  material  of  different  index  in 
one  arm  of  a  Mach-Zehnder  interferometer  will  cause  a  change 
in  phase  velocity  with  respect  to  the  wave  in  the  other  arm. 
The  phase  shift  caused  by  this  change  is 


f  f 

d<s  =  —  -  —  dy 


Cl.  5) 


where  vD  and  vv  are  the  phase  velocities  of  the  waves  in  the 
medium  and  in  the  vacuum,  respectively. 

Equation  Cl. 5)  can  be  simplified  by  noting  that 


Cl.  6) 


C I.  7) 


3y  using  Equations  Cl. 6)  and  Cl. 7)  in  Equation  Cl. 5)  and  then 
integrating  on  both  sides,  the  phase  shift  Cfringe  shift)  due 
to  the  inserted  material  is  found: 


6  =  —  [pCy )- 11  dy  . 

X  JL 


C  I.  8) 
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APPENDIX  II.  APPROXIMATIONS  FOR  PLASMA  INDEX 

The  purpose  of  this  Appendix  is  to  demonstrate  that 
the  refractive  index  for  the  dense  plasma  focus  can  be 
approximated  by 


(ij. 


H  =  (1- 


2cib 


C II.  I) 


which  is  effectively  the  sarne  as  Equation  (2)  in  Chapter  I, 
except  that  the  radial  wave  frequency  is  used  Cu=2jrf).  The 
major  effects  on  the  plasma  index  are  from  the  electron 
collision  frequency  with  heavy  particles  (ve),  the  electron 
cyclotron  frequency  and  the  plasma  frequency  Cup}.  It 
is  also  of  concern  that  the  wave  frequency  is  not  near  the 
upper  hybrid  resonance  frequency,  Cuuh).  Note  that  if  the 
following  approximations  hold  for  the  electrons,  they  will 
also  hold  true  for  the  ions  due  to  their  larger  mass. 

The  electron  collision  frequency  can  be  written^ 


ve  =  2.91’  10-6  neXeiTe_1-5 


(II. 2) 


where  the  Coulomb  logarithm  for  electron-ion  collisions,  Xe^, 
is  given  by 


Xei  =  24-ln(ne*Te-i)  . 
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1 


C II. 3) 


and  an  electron 


For  a  maximum  density3  of  ne=2- 101<3cm  3 


temperature'3 
ve~4‘ 10i0/sec. 


of  500  eV,  the  collision  frequency  is 
The  radial  wave  frequency  (N2)  is 


o=5.  6‘ lO^/sec  and  the  plasma  frequency  for  the  given  density 


is  Up=2.  5  1014/sec. 


The  conditions 


and 


ve 2  <<  (j2Co2-wp2)2/(jp4 


C II. 4) 


C II. 5) 


are  satisfied,  and  the  refractive  index  can  be  approximated 
by6 


2 


(II. 6) 


The  binomial  expansion  can  be  used  to  arrive  at  Equation 
(II. 1),  since  Op2  is  much  less  than  o2  (typically  Op  <  O.lu), 
Since  a  magnetic  field  is  associated  with  the 
pinch,  the  plasma  is  anisotropic.  In  other  words,  the 
refractive  index  is  a  function  of  both  the  direction  of 
propagation  and  electric  field  orientation  of  the  wave 
relative  to  the  magnetic  field  lines.  In  order  to 
demonstrate  that  the  magnetic  field  has  a  negligible  effect 
on  the  index  in  the  case  of  the  plasma  pinch,  the  electron 


cyclotron  frequency  must  first  be  calculated  where 


=  eB/me  •  C II. 7) 

The  magnetic  induction,  B,  can  be  approximated  by  applying 
Ampere's  law  to  an  infinitly  long  wire  with  current  ID.  Then 
the  magnetic  induction  at  the  edge  of  a  pinch  of  radius  R  can 
be  approximated  by 


^o*o 

2nR 


(II.  8) 


The  pinch  current  is  typically  -500  RA.  The  corresponding 
magnetic  induction  for  a  pinch  radius3  of  1  mm  Cmaximum 
compression)  is  B~1  MGauss.  The  conditions  set  an  upper 
limit  for  the  electron  cyclotron  frequency,  0^2' 1013/sec. 

A  wave  propagating  parallel  to  the  field  lines 
becomes  circularly  polarized  with  either  left-  or  right-hand 
rotation.  The  corresponding  index  of  refraction  isb 

9 

M1(  r  =  I  1- - - -  C  II.  9) 

oCo±ob') 

Since  Equation  (II.  1)  is  a  good  approximation. 

If  the  wave  is  propagating  perpendicular  to  the 
magnetic  field  lines,  the  field  has  no  effect  when  the 


electric  vector  of  the  wave  is  parallel  to  it.  The  equation 
in  this  case  is  simply  Equation  (II.  6).  When  the  electric 
vector  is  perpendicular  to  the  magnetic  field  lines 
(extraordinary  wave),  the  index  of  refraction  can  be  written 


p  =  U-  p—  P 
'ex  90  9  9 


(II. 10) 


9  9  o 

and  since  Equation  (II.  1)  is  also  a  good 

approximation  for  Equation  (II. 10). 

Finally,  it  is  noted  that  the  wave  frequency  is 


well  away  from  the  upper  hybrid  frequency  defined  by 


APPENDIX  III.  INTERFEROMETER  ALIGNMENT  PROCEDURE 

The  purpose  of  this  Appendix  is  to  outline  the 
procedure  used  in  alignment  of  the  Mach-Zehnder  inter¬ 
ferometer  of  this  thesis.  Refer  to  Figures  8,  11  and  13  of 
Chapter  II  when  following  the  steps. 
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1.  Align  the  HeNe  laser  beam  through  the  center  of  the 
two  focus  chamber  windows. 

2.  Without  moving  the  HeNe  laser,  align  the  N2  laser 
such  that  the  HeNe  beam  is  concentric  within  its 
window  and  mirror.  This  assures  that  the  two  beams 
will  take  the  same  path. 

3.  Using  the  bolts  that  attach  the  gimbal  mounts  to 
the  focus  flanges,  adjust  BS^  and  M2  to  ~45°  angles 
with  respect  to  the  HeNe  beam  such  that  they  are 
also  centered  about  the  beam. 

4.  The  gimbal  mount  micrometer  adjustment  screws  can 
now  be  used  to  vertically  align  the  HeNe  reference 
and  probe  beams  with  reference  to  a  plumb  bob.  The 
purpose  of  the  plumb  bob  is  to  keep  both  the 
reference  and  probe  beams  in  one  plane  (not  skew). 

5.  Use  the  micrometer  adjustment  screw  of  Mi  to 
shift  the  reference  beam  horizontally  through  BS2 
such  that  it  lines  up  with  the  proi  beam  in  a 
vertical  plane  at  the  output  of  3S-,.  Use  the  other 
micrometer  adjustment  screw  of  to  coextend  the 
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reference  beam  with  the  8"  steel  channel  used  in 
the  detection  stage- 

6.  The  flange  screw  is  now  used  to  position  BS2  to 
%45°  with  respect  to  and  centered  around  the  HeNe 
reference  beam-  The  film  detector  should  now  be 
horizontally  shifted  such  that  HeNeref  is  centered 
within  the  lens  and  bandpass  filter. 

7.  Adjust  Mo  so  that  the  probe  beam  and  the  reference 
beam  coincide  on  the  output  surface  of  BS-?. 

8.  Adjust  BS-,  so  that  “the  beams  coincide  near  the 
pinhole.  Fringes  in  the  HeNe  beam  should  occur. 
Check  near  the  lens. 

9.  Readjust  BS?  to  obtain  horizontal  and  adequately 
spaced  fringes. 


The  next  step  involves  setting  the  path  difference 
between  the  reference  and  probe  waves  to  zero.  Since  the 
nitrogen  laser  used  in  this  experiment  is  a  single  shot  de¬ 
vice,  a  photograph  must  be  usea  after  each  adjustment  of  the 
translational  stage  until  highly  visible  fringes  are  ob¬ 
served.  It  is  suggested  that  the  TS  be  adjusted  in  3  full 
turn  increments  (~1  mm  steps)  until  highly  visible  fringes 
occur  around  the  edges  of  the  fringe  pattern.  Then  the 
increments  should  be  reduced  to  1  turn  until  the  best 
possible  fringes  are  observed.  By  further  reducing  the 
number  of  turns  to  fractional  amounts,  the  quality  of  fringes 
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